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The question of first-cause has troubled philosophers and cosmologists alike. Now that it is 
apparent that our universe began in a Big Bang explosion, the question of what happened before the 
Big Bang arises. Inflation seems like a very promising answer, but as Borde and Vilenkin have shown, 
the inflationary state preceding the Big Bang could not have been infinite in duration — it must have 
had a beginning also. Where did it come from? Ultimately, the difficult question seems to be how to 
make something out of nothing. This paper explores the idea that this is the wrong question — that 
that is not how the Universe got here. Instead, we explore the idea of whether there is anything in 
the laws of physics that would prevent the Universe from creating itself. Because spacetimes can be 
curved and multiply connected, general relativity allows for the possibility of closed timelike curves 
(CTCs). Thus, tracing backwards in time through the original inflationary state we may eventually 
encounter a region of CTCs — giving no first-cause. This region of CTCs may well be over by 
now (being bounded toward the future by a Cauchy horizon). We illustrate that such models — 
with CTCs — are not necessarily inconsistent by demonstrating self-consistent vacuums for Misner 
space and a multiply connected de Sitter space in which the renormalized energy-momentum tensor 
does not diverge as one approaches the Cauchy horizon and solves Einstein's equations. Some 
specific scenarios (out of many possible ones) for this type of model are described. For example: a 
metastable vacuum inflates producing an infinite number of (Big-Bang-type) bubble universes. In 
many of these, either by natural causes or by action of advanced civilizations, a number of bubbles 
of metastable vacuum are created at late times by high energy events. These bubbles will usually 
collapse and form black holes, but occasionally one will tunnel to create an expanding metastable 
vacuum (a baby universe) on the other side of the black hole's Einstein-Rosen bridge as proposed 
by Farhi, Guth, and Guven. One of the expanding metastable-vacuum baby universes produced in 
this way simply turns out to be the original inflating metastable vacuum we began with. We show 
that a Universe with CTCs can be stable against vacuum polarization. And, it can be classically 
stable and self-consistent if and only if the potentials in this Universe are retarded — which gives 
a natural explanation of the arrow of time in our universe. Interestingly, the laws of physics may 
allow the Universe to be its own mother. 

PACS number(s): 98.80.-k 



I. INTRODUCTION 



The question of first-cause has been troubling to philosophers and scientists alike for over two thousand years. 
Aristotle found this sufficiently troubling that he proposed avoiding it by having the Universe exist eternally in both 
the past and future. That way, it was always present and one would not have to ask what caused it to come into being. 
This type of model has been attractive to modern scientists as well. When Einstein developed general relativity and 
applied it to cosmology, his first cosmological model was the Einstein static universe, which had a static S 3 spatial 
geometry which lasted forever, having no beginning and no end jjj. 

As we shall discuss, since the Big Bang model's success, models with a finite beginning have taken precedence, 
even when inflation and quantum tunneling are included. So the problem of first-cause reasserts itself. The big 
question appears to be how to create the universe out of nothing. In this paper we shall explore the idea that this is 
the wrong question. A remarkable p roper ty of general relativity is that it allows solutions that have closed timelike 
curves (CTCs) (for review see |9|,|l0f ). Often, the beginning of the universe, as in Vilenkin's tunneling model 

jllj and Hartle and Hawking's no-boundary model is pictured as being like the south pole of the earth and 
it is usually said that asking what happened before that is like asking what is south of the south pole |13) . But, 
suppose the early universe contains a region of CTCs. Then, asking what was the earliest point might be like asking 
what is the easternmost point on the Earth. You can keep going east around and around the Earth — there is no 
eastern-most point. In such a model every event in the early universe would have events that preceded it. This 
period of CTCs could well have ended by now, being bounded by a Cauchy horizon. Some initial calculations of 
vacuum polarization in spacetimes with CTCs indicated that the renormalized energy-momentum tensor diverged at 
the Cauchy horizon separating the region with CTCs from the region without closed causal curves, or at the polarized 
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hypersurfaces nested inside the Cauchy horizon p4|-p^[ . Some of these results motivated Hawking [ ^9pc| ] to propose 
the chronology protection conjecture which states that the laws of physics do not allow the appearance of CTCs. But, 
a number of people have challenged the chronology protection conjecture by giving counter-examples [p^|j9pl| E$9|. In 
particular, Li and Gott |3(|l have recently found that there is a self-consistent vacuum in Misner space for which the 
renormalized energy-momentum tensor of vacuum polarization is zero everywhere. (Cassidy ]3l| | has independently 
given an existence proof that there should be a quantum state for a conformally coupled scalar field in Misner space, 
for which the renormalized energy-momentum tensor is zero everywhere, but he has not shown what state it should 
be. Li and Gott ]3C|| have found that it is the "adapted" Rindler vacuum.) In this paper we give some examples to 
show how it is possible in principle to find self-consistent vacuum states where the renormalized energy-momentum 
tensor does not blow up as one approaches the Cauchy horizon. To produce such a region of CTCs, the universe 
must, at some later time, be able to reproduce conditions as they were earlier, so that a multiply connected solution 
is possible. Interestingly, inflation is well suited to this. A little piece of inflationary state expands to produce a large 
volume of inflationary state, little pieces of which resemble the starting piece. Also there is the possibility of forming 
baby universes at late times where new pieces of inflating states are formed. Farhi, Guth, and Guven [g2|, Harrison 
f33| , Smolin |34|]35f l , and Garriga and Vilenkin |$(| have considered such models. If one of those later inflating pieces 
simply turns out to be the inflating piece that one started out with, then the Universe can be its own mother. Since 
an infinite number of baby universes are created, as long as the probability of a particular multiple connection forming 
is not exactly zero, then such a connection might be expected, eventually. Then the Universe neither tunneled from 
nothing, nor arose from a singularity; it created itself (Fig. |l|). 

Before discussing this approach to the first-cause problem, let us review just how troublesome this problem has 
been. As we have noted, Einstein Q initially tried to avoid it by siding with Aristotle in proposing a model which 
had an infinite past and future. The Einstein static universe appears to be the geometry Einstein found a priori 
most aesthetically appealing, thus presumably he started with this preferred geometry and substituted it into the 
field equations to determine the energy-momentum tensor required to produce it. He found a source term that 
looks like dust (stars) plus a term that was proportional to the metric which he called the cosmological constant. 
The cosmological constant, because of its homogeneous large negative pressure, exerts a repulsive gravitational effect 
offsetting the attraction of the stars for each other; allowing a static model which could exist (ignoring instabilities, 
which he failed to consider) to the infinite past and future. If one did not require a static model, there would be no need 
for the cosmological constant. Friedmann [p7| c alculated models without it, of positive, negative or zero curvature, 
all of which were dynamical. When Hubble |[38| discovered the expansion of the universe, Einstein pronounced the 
cosmological constant the biggest blunder of his life. 

But now there was a problem: all three Friedmann models (k = 0, k = 1, and k = —I) that were expanding at the 
present epoch had a beginning in the finite past (see e.g. |39| , [40[ | ) . In the Friedmann models the universe began in a 
singularly dense state at a finite time in the past. The equations could not be pushed beyond that finite beginning 
singularity. Furthermore, if today's Hubble constant is Ho, then all of the Friedmann models had ages less than 
t H = Hq 1 . The universe thus began in a Big Bang explosion only a short time ago, a time which could be measured 
in billions of years. The universe was not infinitely old. Gamow [^T]j4^] and his colleagues Alphcr and Herman ^] 
calculated the evolution of such a Big Bang cosmology, concluding correctly that in its early phases it should have 
been very dense and very hot, and that the thermal radiation present in the early universe should still be visible today 
as microwave radiation with a temperature of approximately 5K. Penzias and Wilson's discovery of the radiation 
with a temperature of 2.7K cinched the case for the Big Bang model. The COBE results which have shown 
a beautifully thermal spectrum (4|Q] and small fluctuations in the temperature ST/T = ICU 5 Q, fluctuations 
that are of approximately the right magnitude to grow into the galaxies and clusters of galaxies that we see at the 
present epoch, have served to make the Big Bang model even more certain. With the Big Bang model in ascendancy, 
attention focused on the initial singularity. Hawking and Penrose proved a number of singularity theorems 
showing that, with some reasonable constraints on the energy-momentum tensor, if Einstein's equations are correct 
and the expansion of the universe is as observed today, there is no way to avoid an initial singularity in the model; 
that is, initial singularities would form even in models that were not exactly uniform. So the initial singularity was 
taken to be the first-cause of the Universe. This of course prompted questions of what caused the singularity and 
what happened before the singularity. The standard answer to what happened before the Big Bang singularity has 
been that time was created at the singularity, along with space, and that there was no time before the Big Bang. 
Asking what happened before the Big Bang was considered to be like asking what is south of the south pole. But 
particularly troublesome was the question of what caused the initial singularity to have its almost perfect uniformity 
— for otherwise the microwave background radiation would be of vastly different temperatures in different directions 
on the sky. Yet the initial singularity could not be exactly uniform, for then we would have a perfect Friedmann 
model with no fluctuations which would form no galaxies. It needed to be almost, but not quite perfectly uniform — 
a remarkable situation — how did it get that way? These seemed to be special initial conditions with no explanation 
for how they got that way. 
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Another problem was that singularities in physics are usually smeared by quantum effects. As we extrapolated 
back toward the initial singularity (of infinite density), we would first reach a surface where the density was equal 
to the Planck density and at this epoch classical general relativity would break down. We could not extrapolate 
confidently back to infinite density, we could only say that we would eventually reach a place where quantum effects 
should become important and where classical general relativity no longer applied. Since we do not have a theory of 
quantum gravity or a theory-of-everything we could honestly say that the singularity theorems only told us that we 
would find regions in the early universe where the density exceeded the GUT or Planck densities beyond which we 
did not know what happened — rather much like the Terra Incognita of old maps. We could not then say how our 
universe formed. 

So, questions about how the initial Big Bang singularity was formed and what preceded it remained. The closed 
Friedmann model, popular because it is compact and therefore needs no boundary conditions, re-collapses in a finite 
time in the future to form a Big Crunch singularity at the end. Singularity theorems tell us that in a collapsing 
universe the final Big Crunch singularity cannot be avoided. Classical general relativity tells us that a closed universe 
begins with a singularity and ends with a singularity, with nothing before and nothing after. Nevertheless, many 
people speculated that there could be more than one connected cycle — after all, the singularities only indicated 
a breakdown of classical general relativity and the quantum Terra Incognita at the Planck density might allow a 
cosmology collapsing toward a Big Crunch to bounce and make another Big Bang |5fj|-|5^] . In support of this is the 
fact that de Sitter space (representing the geometry of a false vacuum — an inflationary state as proposed by Guth 
|54f — with a large cosmological constant) looks like a spatially closed S 3 universe whose radius as a function of 
proper time is a(t) = r cosh(t/r ), where r = (3/A) 1 / 2 is the radius of the de Sitter space and A is the cosmological 
constant (throughout the paper we use units G = c = 7j = A;b = 1), which is a collapsing cosmology which bounces 
and turns into an expanding one. Thus if quantum gravitational effects make the geometry look like de Sitter space 
once the density reaches the Planck density as some have suggested ]5S|-|57|], then a Big Crunch singularity might 
be avoided as the closed universe bounced and began a Big Bang all over again. This bouncing model avoids the 
first-cause problem. The answer to what caused our universe in this model is "the collapse of the previous universe" , 
and so on. An infinite number of expansion and contraction cycles make up the Universe (note the capital U — in 
this paper this denotes the ensemble of causally connected universes) which consists of an infinite number of closed 
Big Bang models laid out in time like pearls on a string. The Universe (the infinite string of pearls) has always been 
in existence and will always be in existence, even though our cycle, our standard closed Big Bang cosmology (our 
pearl) has a finite duration. So we are back to Aristotle, with an eternal Universe, and close to Einstein with just an 
oscillating (rather than static) closed Universe that has infinite duration to the past and future. Thus in this picture 
there is no first-cause because the Universe has existed infinitely far back in the past. 

The oscillating universe was thought to have some problems with entropy J53j . Entropy is steadily increasing with 
time, and so each cycle would seem to be more disordered than the one that preceded it. Since our universe has a 
finite entropy per baryon it was argued, there could not be an infinite number of cycles preceding us. Likewise it was 
argued that each cycle of the universe should be larger than the preceding one, so if there were an infinite number 
preceding us, our universe would have to look indistinguishable from flat (i.e., closed but having an infinite radius of 
curvature). The real challenge in this model is to produce initial conditions for our universe (our pearl) that were 
as uniform and low entropy as observed. COBE tells us that our universe at early times was uniform to one part in 
a hundred thousand [ [47|] . At late times we expect the universe at the Big Crunch to be very non-uniform as black 
hole singularities combine to form the Big Crunch. In the early universe the Weyl tensor is zero, whereas at the Big 
Crunch it would be large ]5S| , |59| ]. How does the chaotic high-entropy state at the Big Crunch get recycled into the 
low-entropy, nearly uniform, state of the next Big Bang? If it does not, then after an infinite number of cycles, why 
are we not in a universe with chaotic initial conditions? 

Entropy and the direction of time may be intimately tied up with this difference between the Big Bang and the 
Big Crunch. Maxwell's equations (and the field equations of general relativity) are time-symmetric, so why do we 
see only retarded potentials? Wheeler and Feynman addressed this with their absorber theory [|BCj . They supposed 
that an electron shaken today produces half-advanced-half-retarded fields. The half-advanced fields propagate back 
in time toward the early universe where they are absorbed (towards the past the universe is a perfect absorber) 
by shaking charged particles in the early universe. These charged particles in turn emit half-advanced-half-retarded 
fields; their half-retarded fields propagate toward the future where they: (a) perfectly cancel the half-advanced fields 
of the original electron, (b) add to its retarded fields to produce the electron's full retarded field, and (c) produce a 
force on the electron which is equal to the classical radiative reaction force. Thus, the electron only experiences forces 
due to fields from other charged particles. This is a particularly ingenious solution. It requires only that the early 
universe is opaque — which it is — and that the initial conditions are low-entropy; that is, there is a cancelation of 
half-advanced fields from the future by half-retarded fields from the past, leaving no "signals" in the early universe 
from later events — a state of low-entropy. (Note that this argument works equally well in an open universe where the 
universe may not be optically thick toward the future — all that is required is that the universe be a perfect absorber 
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in the past, i.e., toward the state of low-entropy.) Wheeler and Feynman noted that entropy is time-symmetric like 
Maxwell's equations. If you find an ice cube on the stove, and then come back and re-observe it a minute later, you 
will likely find it half- melted. Usually an ice cube gets on a stove by someone just putting it there (initial conditions), 
but suppose we had a truly isolated system so that the ice cube we found was just a statistical fluctuation. Then if 
we asked what we would see if we had observed one minute before our first observation, we will also be likely to see a 
half-melted ice cube, for finding a still larger ice cube one minute before would be unlikely because it would represent 
an even more unlikely statistical fluctuation than the original ice cube. In an isolated system, an (improbable) state of 
low-entropy is likely to be both followed and preceded by states of higher-entropy in a time-symmetric fashion. Given 
that the early universe represents a state of hig h order, it is thus not surprising to find entropy increasing after that. 
Thus, according to Wheeler and Feynman |S0|, the fact that the retarded potentials arrow of time and the entropy 
arrow of time point in the same direction is simply a reflection of the low-entropy nature of the Big Bang. The Big 
Crunch is high-entropy, so time follows from past to future between the Big Bang and the Big Crunch. 

Thus, in an oscillating universe scenario, we might expect entropy to go in the opposite direction with respect to 
time, in the previous cycle of oscillation. In that previous universe there would be only advanced potentials and 
observers there would sense a direction of time opposite to ours (and would have a reversed definition of matter and 
anti- matter because of CPT invariance). Thus the cycle previous to us would, according to our definition of time, have 
advanced potentials and would end with a uniform low-entropy Big Crunch and begin with a chaotic high-entropy 
Big Bang (see Gott |Sl| for further discussion). Thus, an infinite string of oscillating universes could have alternating 
high and low-entropy singularities, with the direction of the entropy (and causality — via electromagnetic potentials) 
time-reversing on each succeeding cycle. Every observer using the entropy direction of time would see in his "past" 
a low-entropy singularity (which he would call a Big Bang) and in his "future" a high-entropy singularity (which he 
could call a Big Crunch). Then the mystery is why the low-entropy Big Bangs exist — they now look improbable. 
An oscillating universe with chaotic bangs and crunches and half-advanced-half-retarded potentials throughout would 
seem more likely. At this point anthropic arguments p3] could be brought in to say that only low-entropy Big Bangs 
might produce intelligent observers and that, with an infinite number of universes in the string, eventually there 
would be — by chance — a sufficiently low-entropy Big Bang to produce intelligent observers. Still, the uniformity 
of the early universe that we observe seems to be more than that required to produce intelligent observers, so we 
might wonder whether a random intelligent observer in such a Universe would be expected to see initial conditions in 
his/her Big Bang as uniform as ours. (Among intelligent observers, the Copernican principle tells us that you should 
not expect to be special. Out of all the places for intelligent observers to be there are by definition only a few special 
places and many non-special places, so you should expect to be in one of the many non-special places 



II. INFLATION AS A SOLUTION 



Guth's proposal of inflation [|54| offered an explanation of why the initial conditions in the Big Bang should be 
approximately, but not exactly uniform. (For review of inflation see p4] , p5| .) In the standard Big Bang cosmology 
this was always a puzzle because antipodal points on the sky on the last scattering surface at 1 + z ~ 1000 had not 
had time to be in communication with each other. When we see two regions which are at the same temperature, 
the usual explanation is that they have at some time in the past been in causal communication and have reached 
thermal equilibrium with each other. But there is not enough time to do this in the standard Big Bang model where 
the expansion of the scale factor at early times is a(t) cx t x l 2 . Grand unified theories (GUT) of particle physics 
suggest that at early times there might have been a non-zero cosmological constant A, which then decayed to the zero 
cosmological constant we see today. This means that the early universe approximates de Sitter space with a radius 
r = (3/A) 1 / 2 whose expansion rate at late times approaches a(t) = roexp(i/Vo). Regions that start off very close 
together, and have time to thermally equilibrate, end up very far apart. When they become separated by a distance 
ro, they effectively pass out of causal contact — if inflation were to continue forever, they would be beyond each 
other's event horizons. But eventually the epoch of inflation ends, the energy density of the cosmological constant 
is dumped into thermal radiation, and the expansion then continues as a(t) cx t 1 / 2 as in a radiation-dominated Big 
Bang cosmology. As the regions slow their expansion from each other, enough time elapses so that they are able to 
interchange photons once again and they come back into effective causal contact. As Bill Press once said, they say 
"hello" , "goodbye" , and "hello again" . When they say "hello again" they appear just like regions in a standard Big 
Bang cosmology that are saying "hello" for the first time (i.e., are just coming within the particle horizon) except 
that with inflation these regions are already in thermal equilibrium with each other, because they have seen each 
other in the past. Inflation also gives a natural explanation for why the observed radius of curvature of the universe 
is so large (a > cHq 1 ~ 3000/i _1 Mpc; here Hq = 100ft. km s _1 Mpc -1 is the Hubble constant). During the Big 
Bang phase, as the universe expands, the radius of the universe a expands by the same factor as the characteristic 
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wavelength A of the microwave background photons, so a/ A = costant > e 67 . How should we explain this large 
observed dimensionless number? Inflation makes this easy. The energy density during the inflationary epoch is A/8tt. 
Let A be the characteristic wavelength of thermal radiation which would have that density. Even if a started out of 
the same order as A, by the end of the inflationary epoch a > Ae 67 , providing that the inflationary epoch lasts at least 
as long as 67ro, or 67 e-folding times. At the end of the inflationary epoch when the inflationary vacuum of density 
A/8-7T decays and is converted into an equivalent amount of thermal radiation, the wavelength of that radiation will 
be A and the ratio of a/ A is fixed at a constant value which is a dimensionless constant > e 67 , retained as the universe 
continues to expand in the radiation and matter-dominated epochs. Thus, even a short run of inflation, of 67 e-folding 
times or more, is sufficient to explain why the universe is as large as it is observed to be. 

Another success of inflation is that the observed Zeldovich-Peebles-Yu-Harrison fluctuation spectrum with index 
n = 1 pq-plfl has been naturally predicted as the result of random quantum fluctuations |]69|-[72[. The inflationary 
power spectrum with CDM has been amazingly successful in explaining the qualitative features of observed galaxy 
clustering (cf. J73|-p2| ) . The amount of large scale power seen in the observations suggests an inflationary CDM power 
spectrum with 0.2 < Qh < 0.3 f§|-§|]. 



III. OPEN BUBBLE UNIVERSES 

Gott fi^l has shown how an open inflationary model might be produced. The initial inflationary state approximates 
de Sitter space, which can be pictured by embedding it as the surface W 2 + X 2 +Y 2 + Z 2 — V 2 — r\ in a five-dimensional 
Minkowski space with metric ds 2 = -dV 2 + dW 2 + dX 2 + dY 2 + dZ 2 Slice de Sitter space along surfaces of 

V = constant, then the slices are three-spheres of positive curvature W 2 + X^ + Y 2 + Z 2 = a 2 where a 2 = r 2 . + V 2 . If 
t measures the proper time, then V = ro sinh(i/ro) and a(t) — ro cosh(£/ro). This is a closed universe that contracts 
then re-expands — at late times expanding exponentially as a function of proper time. If slices of V + X = constant 
are chosen, the slices have a flat geometry and the expansion is exponential with a(t) = ro exp(t/ro). If the slices are 
vertical (W = constant > ro), then the intersection with the surface is if 3 , a hyperboloid X 2 + Y 2 + Z 2 — V 2 = —a 2 
living in a Minkowski space, where a 2 — W 2 — r 2 . This is a negatively curved surface with a radius of curvature a. 
Let t be the proper time from the event E (W = ro, X — 0, Y = 0, Z = 0, V = 0) in the de Sitter space. Then the 
entire future of E can be described as an open k = —1 cosmology where a(t) — ro sinh(i/ro). At early times, t <C ro, 
near E, a(t) oc t, and the model resembles a Milne cosmology [£1[, but at late times the model expands exponentially 
with time as expected for inflation. This is a negatively curved (open) Friedmann model with a cosmological constant 
and nothing else. Note that the entire negatively curved hyperboloid (H 3 ), which extends to infinity, is nevertheless 
causally connected because all points on it have the event E in their past light cone. Thus, the universe should have 
a microwave background that is isotropic, except for small quantum fluctuations. At a proper time t\ after the event 
E, the cosmological constant would decay leaving us with a hot Big Bang open (fc = —I) cosmology with a radius of 
curvature of a = ro sinh(ri/ro) at the end of the inflationary epoch. If t\ = 67ro, then f2 is a few tenths today; if 
ti > 67r , then ~ I today 

Gott jS 7 ^ noted that this solution looks just like the interior of a Coleman bubble p2|| . Coleman and de Luccia |93| ] 
showed that if a metastable symmetric vacuum (with the Higgs field = 0), with positive cosmological constant A 
were to decay by tunneling directly through a barrier to reach the current vacuum with a zero cosmological constant 
(where the Higgs field ip — <f>o), then it would do this by forming a bubble of low-density vacuum of radius a around 
an event E. The pressure inside the bubble is zero while the pressure outside is negative (equal to — A/87r), so the 
bubble wall accelerates outward, forming in spacetime a hyperboloid of one sheet (a slice of de Sitter space with 
W = constant < ro). This bubble wall surrounds and is asymptotic to the future light cone of E. If the tunneling is 
direct, the space inside the bubble is Minkowski space (like a slice W = constant < ro in the embedding space, which 
is flat). The inside of the future light cone of E thus looks like a Milne cosmology with Q = and a(t) = t. Gott |89| ] 
noted that what was needed to produce a realistic open model with O of a few tenths today was to have the inflation 
continue inside the bubble for about 67 e-folding times. Thus, our universe was one of the bubbles and this solved 
the problem of Guth's inflation that in general one expected the bubbles not to percolate |9lJ9j|. But, from inside 
one of the bubbles, our view could be isotropic Q. 

It was not long before a concrete mechanism to produce such continued inflation inside the bubble was proposed. 
A couple of weeks after Gott's paper appeared Linde's |9(| proposal of new inflation appeared, followed shortly by 
Albrecht and Steinhardt They proposed that the Higgs vacuum potential V((j>) had a local minimum at <fi — 
where 1^(0) = A/8ir. Then there was a barrier at (j> = (j>±, followed by a long flat plateau from (f>\ to c/)q where it drops 
precipitately to zero at (ftp. The relation of this to the open bubble universe's geometry is outlined by Gott (9F|] (see 
Fig. I and Fig. 2 in l^]). The de Sitter space outside the bubble wall has 4> = 0. Between the bubble wall, at a 
spacelike separation a from the event E, and the end of the inflation at the hyperboloid H 3 which is the set of points 
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at a future timelike separation of t\ from E, the Higgs field is between <pi and </>o, and t\ is the time it takes the 
field (after tunneling) to roll along the long plateau [where V((f>) is approximately equal to A/8-7T and the geometry is 
approximately de Sitter] . After that epoch, <fi = where the energy density has been dumped into thermal radiation 
and the vacuum density is zero (i.e., a standard open Big Bang model). In order that inflation proceeds and the 
bubbles do not percolate, it is required that the probability of forming a bubble in de Sitter space per four volume r$ 
is e < e cr where 5.8 x 1CP 9 < e cr < 0.24 fl95| . In order that there be a greater than 5% chance that no bubble should 
have collided with our bubble by now, so as to be visible in our past light cone, e < 0.01 for SI = 0.4, A = 0, h = 0.63 
today 1 88 1 , but this is no problem since we expect tunneling probabilities through a barrier to be exponentially small. 
This model has an event horizon, which is the future light cone of an event E' (W = — ro, X = 0, Y — 0, Z = 0, V = 0) 
which is antipodal to E. Light from events within the future light cone of E' never reaches events inside the future 
light cone of E. So we are surrounded by an event horizon. T his produces Hawking radiation; and, if r is of order 
the Planck length, then the Gibbons-Hawking thermal state [ 100 1 (which looks like a cosmological constant due to 
the trace anomaly [101]) should be dynamically important p9| . 

If we observe SI < 1 and SI a = 0, then k = — 1 and we need inflation more than ever — we still need it to explain 
the isotropy of the microwave background radiation and we would now have a large but finite radius of curvature to 
explain, which 67 e-folds of inflation could naturally produce. When Gott told this to Linde in 1982, Linde said, yes, 
if we found that SI < 1, he would still have to believe in inflation but he would have a headache in the morning! Why? 
Because one has to produce a particular amount of inflation, approximately 67 e-folds. If there were 670 e-folds or 670 
million e-folds, then S7 currently would be only slightly less than 1. So there would be what is called a "fine tuning 
of parameters" needed to produce the observed results. 

The single-bubble open inflationary model |8j| discussed above has recently come back into fashion because of a 
number of important developments. On the theoretical side, Ratra and Peebles [102 ,103] have shown how to calculate 
quantum fluctuations in the H 3 hyperbolic geometry with a(t) = rosinh(i/ro) during the inflationary epoch inside 
the bubble in the single bubble m odel. This allows predictions of fluctuations in the microwave background. Bu chcr , 
Goldhaber, and Turok [104, 105 1 have exte nded these calculations, as well as Yamamoto, Sasaki and Tanaka [ 106 1 . 
Importantly, they have explained | 104,|l05| tha t the fine tuning in these models is only "logarithmic" and, therefore, 



not so serious. Linde and Mezhlmian [ 107j 



have shown how there are reasonable potentials whic h cou ld produce 
such bubble universes with different values of SI. In a standard chaotic inflationary potential V(<j>) [109], one could 
simply build in a bump, so that one would randomly walk to the top of the curve via quantum fluctuations and 
then roll down till one lodged behind the bump in a metastable local minimum. One would then tunnel through 
the bump, forming bubbles that would roll down to the bottom in a time T\. One could have a two-dimensional 
potential V(<ft, a) = ^g 2 <p 2 cr 2 + V(a), where g is a constant and there is a metastable trough at a = with altitude 
V(<f>, 0) = A/87T with a barrier on both sides, but one coul d tunnel through the barrier to reach a > where V(<p, a) 
has a true minimum, and at fixed c, is proportional to <f> 2 [107,10?]. Then individual bubbles could tunnel across the 
barrier at different values of 4>, and hence have different roll-down times t\ and thus different values of SI. With a 
myriad of open universes being created, anthropic arguments |6^] come into play and if shorter roll-down times were 
more probable than large ones, we might not be surprised to find ourselves in a model which had SI of a few tenths, 
since if SI is too small, no galaxies will form [ 1 10 1 . 

A second reason for the renaissance of these open inflationary models is the observational data. A number of recent 
estim ates of h (the pre sent Hubble constant i n un its of 100 km s _1 Mpc " 1 ) h ave been made (i.e., h = 0.65 ± 0.06 
| mil , 0.68 < h < 0.77 (m), 0.55 < h < 0.6 1 |lf] , and h = 0.64 ± 0.06 [£]J). Ages of globular cluster stars have 
a 2cr lower limit of about 11.6 billion years [115], we require h < 0.56 if SI = 1, but a more acceptable h < 0.65 if 
SI = 0.4, S7a = 0. Models with low SI but SI + S1a = 1 are also acceptable. Also, studies of large scale structure have 
shown that with the inflationary CDM power spectrum, the standard Q = 1, h = 0.5 model simply does not have 
enough power at large scales. A variety of observational samples and methods have suggested this: counts in cells, 
angular covariance function on the sky, power spectrum analysis of 3D samples, and finally topological analysis, all 
showing that 0.2 < Slh < 0.3 |82-88|. If ft- > 0.55 this implies SI < 0.55, which also agrees with wha t one would 
deduce from the age argument as well as the measured masses in groups and clusters of galaxies 116]. With the 



COBE normalization there is also the problem that with fi = 1, (SM/M) sh - 



L Mpc 



1.1 — 1.5 and this would require 



galaxies to be anti-biased [since for galaxi es {5 M/M ) 8fe -i MDC = 1] and would also lead to an excess of large-separation 
gravitational lenses over those observed 1 1 1 7 ] . These things have forced even enth usias ts of k = models to move 
to models with SI < 1 and a cosmological constant so that SI + SI a = 1 and k — 1 1 1 8 1 . They then have to explain 
the small ratio of the cosmological constant to the Planck density (10~ 120 ). Currently we do not have such a natural 
explanation for a small yet finite A as inflation naturally provides for explaining why the radius of curvature should 
be a big n umb er in the k = — 1 case. 

Turner [119] and Fukugita, Futamase, and Kasai |120| showed that a flat S1a = 1 model produces about 10 times 
as many gravitational lenses as a flat model with SI = 1, and Kochanek [ 1 2 1 1 was able to set a 95% confidence lower 
limit of 0.34 < SI in flat models where SI + S1a = L and a 90% confidence lower limit 0.15 < SI in open models with 
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f^A = 0. Thus, extreme-A dominated models are ruled out by producing too many gravitational lenses. 

Data on cosmic microwave background fluctuations for spherical harmonic modes from I = 2 to I = 500 will provide 
a strong test of these models. With fish 2 = 0.0125, the f2 = 1, f^A = model powe r sp ectrum reaches its peak value 
at I = 200; an £1 = 0.3, Q\ = 0.7 mo del r eaches its peak value also at I = 200 [122|; while an = 0.4, = 
model reaches its peak value at I — 350 H123H . This should be decided by the MAP and PLANCK satellites which will 
measure this range with high accuracy |124[] . 

For the rest of this paper we shall usually assume single-bubble open inflationary models for our Big Bang universe 
(while recognizing that chaotic inflationary models and models with multiple epochs of inflation are also possible; it 
is interesting to note that Penr ose also prefers an open universe from the point of view of the complex- holomorphic 
ideology of his twister theory |l25| ). If the inflation within the bubble is of order 67 e-folds, then we can have Q 
of a few tenths; but if it is longer than that, we will usually see fi near 1 today. In any case, we will be assuming 
an initial metastablc vacuum which decays by forming bubbles through barrier penetration. The bubble formation 
rate per unit four volume r ^ is thus expected to be exponentially small so the bubbles do not percolate. Inflation is 
thus eternal to the future [ 12S| -131]. Borde and Vilenkin have proved that if the Universe were infinitely old (i.e., 
if the de Sit ter spac e were complete) then the bubbles would percolate immediately and inflation would never get 
started (see [ 126 , 127 1 and references cited therein). Recall that a complete de Sitter space may be covered with an S 3 
coordinate system (a k = 1 cosmology) whose radius varies as a(t) = ro cosh(i/ro) so that for early times (t < 0) the 
universe would be contracting and bubbles would quickly collide preventing the inflation from ever reaching t = 0. 
Thus Borde and Vilenkin have proved that in the inflationary scenario the universe must have a beginning. If it starts 
with a three-sphere of radius ro at time t — 0, and after that expands like a{i) — tq cosh(£/ro), the bubbles do not 
percolate (given that the bubble formation rate per four volume Tq is e <C 1) and the inflation continues eternally 
to t = oo producing an infinite number of open bubble universes. Since the number of bubbles forming increases 
exponentially with time without limit, our universe is expected to form at a finite but arbitrarily large time after 
the beginning of the inflationary state. In this picture our universe (our bubble) is only 12 billion years old, but the 
Universe as a whole (the entire bubble forming inflationary state) is of a finite but arbitrarily old age. 



IV. VILENKIN'S TUNNELING UNIVERSE AND HARTLE-HAWKING'S NO-BOUNDARY PROPOSAL 

But how to produce that initial spherical S 3 universe? Vilenkin [TT| suggested that it could be formed from quantum 
tunneling. Consider the embedding diagram for de Sitter space. De Sitter space can be embedded as the surface 
W 2 +X 2 +Y 2 + Z 2 -V 2 = r 2 in a five-dimensional Minkowski space with metric ds 2 = -dV 2 + dW 2 + dX 2 +dY 2 +dZ 2 . 
This can be seen as an S* 3 cosmology with radius a(t) = Tq cosh(i/ro) where V — ro sinh(t/ro) and a 2 — W 2 + X 2 + 
Y 2 + Z 2 gives the geometry of S 3 . This solution represents a classical trajectory with a turning point at a = ro. But 
just as it reaches this turning point it could tunnel to a = where the trajectory may be shown as a hemisphere 
of the Euclidean four-sphere W 2 + X 2 + Y 2 + Z 2 + V 2 = r^ embedded in a flat Euclidean space with the metric 
ds 2 = dV 2 + dW 2 + dX 2 + dY 2 + dZ 2 and a{t E ) = r cos(t E /r ) where a 2 = W 2 + X 2 + Y 2 + Z 2 and V = r sm(t E /r ). 
The time-reversed version of this process would show tunneling from a point at (V = — ro, W = 0, X = 0, Y = 0, Z = 0) 
to a three sphere at V — of radius ro which then expands with proper time like a(t) = ro cosh(t/ro) giving a normal 
de Sitter space — thus Vilenkin's universe created from nothing is obtained p"l"| . 

Hawking has noted that in this case, in Hartle and Hawking's formulation, the point (V = — ro, W = 0, X = 0, Y = 
0, Z = 0) is not special, the curvature does not blow up there: it is like other points in the Euclidean hemispherical 
section [jl3[ . However, this point is still the earliest point in Euclidean time since it is at the center of the hemisphere 
specified by the Euclidean boundary at V — 0. So the beginning point in the Vilenkin model is indeed like the south 
pole of the Earth Q. 

Vilenkin's tunneling universe was based on an analogy between quantum creation of universes and tunneling in 
ordinary quantum mechanics p| . In ordinary quantum mechanics, a particle bounded in a well surrounded by a 
barrier has a finite probabilit y to tunne l through the barrier to the outside if the height of the barrier is finite (as in the 



a-decay of radioactive nuclei [ 132 - 134 1 ) . The wave function outside the barrier is an outgoing wave, the wave function 
in the well is the superposition of an outgoing wave and an ingoing wave which is the reflection of the outgoing wave 
by the barrier. Due to the conservation of current, there is a net outgoing current in the well. The probability for 
the particle staying in the well is much greater than the probability for the particle running out of the barrier. The 
energy of the particle in the well cannot be zero, otherwise the uncertainty principle is violated. Thus there is always 
a finite zero-point-energy. The Vilenkin universe was supposed to be created from "nothing" , where according to 



Vilenkin "nothing" means "a state with no classical spacetime" 135 1. Thus this is essentially different from tunneling 
in ordinary quantum mechanics since in ordinary quantum mechanics tunneling always takes place from one classically 
allowed region to another classically allowed region where the current and the probability are conserved. But creation 
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from "nothing" is supposed to take place from a classically forbidden (Euclidean) region to a classically allowed 
(Lorentzian) region, so the conservation of current is obviously violated. Vilenkin obta ined his tunneling universe by 
choosing a so-called "tunneling boundary condition" for the Wheeler-DeWitt equation 135. 136]. His "tunneling from 
nothing" boundary condition demands that when the universe is big (a 2 A/3 > 1 where A i s the cos mological constant 
and a is the scale factor of the universe) there is only an outgoing wave in the superspace 1 135 136]. If the prob ability 
and current are conserved (in fact there does exist a conserved current for the Wheeler-DeWitt equation [137], and a 
classically allowed solution with a = and zero "energy"), there must be a finite probability for the universe being 
in the state before tunneling (i.e., a = 0) and this probability is much bigger than the probability for tunneling. This 
implies that there must be "something" instead of "nothing" before tunneling. This becomes more clear if matter 
fields are included in considering the creation of universes. In the case of a cosmological constant A and a conformally 
coupled scalar field </> (conformal fields are interesting not only for their simplicity but also because electromagnetic 
fields are conformally invariant) as the source terms in Einstein's equations, in the mini-superspace model (where 
the configurations are the scale factor a of the S 3 Roberts on- W alker metric and a homogeneous conformally coupled 
scalar field </>) the Wheeler-DeWitt equation separates |i~2,138] 



^-^ +X '"') <T>!Xi=Et ™ X! - 
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where \&(a)$(x) is the wave function of the universe [X = (4 7 r/3) 1 / 2 H, e is the "energy level" of the conformally 
coupled scalar field , (we use quotes because for radiation the conserved quantity is E = 4irpa 4 /3 instead of the energy 
47rpa 3 /3 where p is the energy density), and p is a constant determining the operator ordering. Eq. (|l]) is just the 
Schrodinger equation of a harmonic oscillator with unit mass and unit frequency and energy E, the eigenvalues of E 
are n + \ where n = 0, 1, 2, ... Eq. (g) is equivalent to the Schrodinger equation for a unit mass particle with total 
energy E = n + | in the one-dimensional potential 



U(a) = -\ <r 



A 



(3) 



It is clear that in the case of n < — 1), there exist one classically forbidden region a\ < a < a 2 and two classically 

(Fig. H). Because U(a) is regular 



allowed regions < a < a\ and a > a 2 where a{ 2 = 1 "F y 1 — + 
at a — 0, we expect that the wave function 'J (a) is also regular at a = 0. If A <C 1 and the conformally coupled 
scalar field is in the ground state with n = 0, we have a\ ~ 1, a 2 — (3/A) 1 / 2 and the potential in region < a < a\ 
is U{a) ~ ia 2 like a harmonic oscillator. The quantum behavior of the universe in region < a < a\ is like a 
quantum harmonic oscillator. This may describe a quantum oscillating (Lorentzian) universe without Big Bang or 
Big Crunch singularities, which has a finite (but small) probability [~ exp(— 1/A)] to tunnel through the barrier to 
form a de Sitter-type inflating universe. The existence of this tiny oscillating universe is due to the existence of a 
finite "zero-point-energy" (1/2) of a conformally coupled scalar field and this "zero-point-energy" is required by the 
uncertainty principle. Since a conformally coupled scalar field has an equation of state like that of radiation, the 
Friedmann equation for k — +1 is 



C 

7^2 



A 



(4) 



where C = 8tt pa 4 / '3 = constant and p is the energy density of the conformally coupled scalar field. Eq. (||) is equivalent 
to the energy-conservation equation for a classical unit mass particle with zero total energy moving in the potential 



If. A 2 C 



V{a) = - l-_ a '-_ 



(5) 



The difference between U (a) and Via) is caused by the fact that in the integral of action the volume element contains 
a factor a 3 which is also varied when one makes the variation to obtain the dynamical equations. The potential V(a) 
is singular at a = and near a — we have V(a) ~ "20^- F° r A <C 1 and n = (we take C = 2E = 2n + 1), the 
classical universe in region < a < a\ is radiation dominated. This universe expands from a Big Bang singularity, 
reaches a maximum radius, then re-collapses to a Big Crunch singularity: a — is a singularity in the classical picture. 
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But from the above discussion, the Wheeler-DeWitt equation gives a regular wave function at a = 0. In such a case 
near a = the quantum behavior of the universe is different from classical behavior. This implies that, near a = 0, 
classical general relativity breaks down and quantum gravity may remove singularities. This case is like that of a 
hydrogen atom where the classical instability (according to classical electrodynamics, an electron around a hydrogen 
nucleus will fall into the nucleus due to electromagnetic radiation) is cured by quantum mechanics. Anyway, it is 
not nothing at a = 0. There is a small classically allowed, oscillating, radiation dominated, closed, quantum (by 
"quantum" we mean that its quantum behavior deviates significantly from its classical behavior) Friedmann universe 
near a = 0, which has a small probability to tunnel through the barrier to form an inflationary universe. (If A > 0.75 
there is no classically forbidden region and thus no tunneling.) 

So in this model the universe did not come from a point (nothing) but from a tiny classically allowed, oscillating, 
quantum Friedmann universe whose radius is of order the Planck magnitude. But where did this oscillating universe 
come from? Because it has a finite probability to tunnel (each time it reaches maximum radius) to a de Sitter space, 
it has a finite "half-life" for decay into the de Sitter phase and cannot last forever. It could, of course, originate by 
tunneling from a collapsing de Sitter phase (the time-reversed version of the creation of a de Sitter state from the 
oscillating state), but then we are back where we started. In fact, starting with a collapsing de Sitter phase one 
is more likely to obtain an expanding de Sitter phase by simply re-expanding at the classical turning point rather 
than tunneling into and then out of the tiny oscillating universe state. An alternative might be to have the original 
tiny oscillating universe created via a qu antu m fluctuation (since it has just the "zero-point-energy" ) but here we are 



basically returning to the idea of Tryon [139] that you could get an entire Friedmann universe of any size directly via 
quantum fluctuation. But quantum fluctuation of what? You have to have laws of physics and a potential etc. 

Hartle and Hawking [T^ ] made their no-boundary proposal and obtained a model of the universe similar to Vilenkin's 
tunneling universe. The no-boundary proposal is expressed in terms of a Euclidean path integral of the wave function 
of the universe 

9{hab,<f> 1 ,dM)=^2 ( Vg ab V(Pexp[-I(g ab ,cl>,M)}, (6) 

M J 

where the summation is over compact manifolds M with the prescribed boundary dM (being a compact three- 
manifold representing the shape of the universe at a given epoch) as the only boundary; g a b is the Euclidean metric 
on the manifold M with induced three-metric h a b on dM, <fi is the matter field with induced value <fii on dM; I is the 
Euclidean action obtained from the Lorentzian action S via Wick rotation: / = —iS(t — > — ir). In the mini-superspace 
model the configuration space is taken to include the k = +1 Robertson- Walker metric and a homogeneous matter 
field. In the WKB approximation the wave function is (up to a normalization factor) 

* ~ B m exp[-/ cl ( 5afc , 0, M)], (7) 

M 

where I c \ is the Euclidean action for the solutions of the Euclidean field equations (Einstein's equations and matter 
field equations). The factor Bm is the determinant of small fluctuations around solutions of the field equations Jl2| . 
If the matter field is a conformally coupled scalar field <f> = (3/47r) 1 / 2 x/a (which is the case that Hartle and Hawking 
PH discussed), pa 4 is conserved where p is the energy density of (f> satisfying the field equations. Then the Friedmann 
equation is given by Eq. (^). The corresponding Euclidean equation is obtained from Eq. (^) via t — > —it 

da\ 2 A , C 



The solution to Eq. (g) is (for the case |AC < I) 



a(r)=H- 1 



1 -+ 1 -(l-AH 2 C) 1/2 co S (2Hr) 



1/2 



(9) 



1 /2 

where H = (j) 1 ^ 2 - This is a Euclidean bouncing space with a maximum radius a max = H^ 1 [5 + |(1 — 4i7 2 C) 1 / 2 ] 

and a minimum radius et m i n = H^ 1 [5 — |(1 — 4i/ 2 C) 1 / 2 ] (Fig. ||). If C = 0, we have a max = H^ 1 , a m i n = 0, 
and o(t) = H^ 1 cos(Ht), one copy of this bouncing space is a four-sphere with the Euclidean de Sitter metric 
ds 2 = dr 2 +H~ 2 cos 2 (Ht) [d\ 2 + sin 2 x(d6 2 +sin 2 9d(f> 2 )] — which is just a four-sphere embedded in a five-dimensional 
Euclidean space (V, W, X, Y, Z) with metric ds 2 = dV 2 + dW 2 + dX 2 + dY 2 + dZ 2 — this is the solution that Hartle 



and Hawking used p2[. But, as we have argued above, according to Hartle and Hawking |12j and Hawking [138|, the 
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Wheeler-DeWitt equation for ^(x) [Eq. ([[])] gives rise to a "zero-point-energy" for the conformally coupled scalar field: 
Co = 2E(n = 0) = 1 (the state with C — violates the uncertainty principle). One copy of this bouncing Euclidean 
space is not a compact four-dimensional manifold with no boundaries, but has two boundaries with a = a m i n (see 
Fig. |). If H <C 1 (i.e. A < 1), we have a max ~ if -1 , a min ^ 1. 



Penrose [125] has criticized Hawking's no-boundary proposal and the model obtained by gluing a de Sitter space 
onto a four-sphere hemisphere by pointing out that there are only very few spaces for which one can glue a Euclidean 
and a Lorentzian solution together since it is required that they have both a Euclidean and a Lorentzian solution, 
but the generic case is certainly very far from that. Here "with a zero-point-energy" we have have both a Euclidean 
solution and a Lorentzian solution, and they can be glued together. But the Euclidean solution is not closed in 
any way; that is, it does not enforce the no-boundary proposal. Hartle and Hawking argued that there should be a 
constant eo m E which arises from the renormalization of the matter field, i.e., E should be n+i + £ p2| . But there 
is no reason that eo should be — ^ to exactly cancel the "zero-point-energy" ^ . (As in the case of a quantum harmonic 
oscillator, we have no reason to neglect the zero-point-energy.) In fact, since eo comes from the renormalization of 
the matter field (without quantization of gravity), it shou ld be much less than the Planck magnitude, i.e., eo 1, 



and thus eo is negligible compared with i. In fact in [138 Hawking has dropped eo. 

In |1^] Hartle and Hawking have realized that for excited states (n > 0), there are two kinds of classical solutions: 
one represents universes which expand from zero volume, to reach a maximum radius, and then re-collapse (like our 
tiny oscillating universe); the other represents the de Sitter- type state of continual expansion. There are probabilities 
for a universe to tunnel from one state to the other. Here we argue that for the ground state (n = 0), there are 
also two such kinds of Lorentzian universes. One is a tiny quantum oscillating universe (having a maximum radius 
with Planck magnitude). Here "quantum" just means that the classical description fails (so singularities might be 
removed). The other is a big de Sitter- type universe. These two universes can be joined to one another through a 
Euclidean section, which describes quantum tunneling from a tiny oscillating universe to an inflating universe (or from 
a contracting de Sitter- type universe to a tiny oscillating universe). During the tunneling, the radius of the universe 
makes a jump (from the Planck length to H^ 1 or vice versa). 

As Hartle and Hawking |l^] calculated the wave function of the universe for the ground state, they argued that, for 
the conformally coupled scalar field case, the path integral over a and \ = (47r/3) 1 / 2 0a separates since "not only the 
action separates into a sum of a gravitational part and a matter part, but the boundary condition on the a(j]) and 
x(r?) summed over do not depend on one another" where r\ is the conformal time. The critical point for the variable's 
separation in the path integral is that "the ground state boundary conditions imply that geometries in the sum are 
conformal to half of a Euclidean-Einstein static universe; i.e., the range of r\ is (— oo,0). The boundary conditions 
at infinite r\ are that x( 7 ?) an d o,(rj) vanish. The boundary conditions at 77 = are that a(0) and x(0) match the 
arguments of the wave function ao and xo" flj]- But this holds only for some specific cases, such as de Sitter space. 
Our solution (j^j does not obey Hartle and Hawking's assumption that r] ranges from —00 to 0. For a general k = +1 
(Euclidean) Robertson- Walker metric, rj — J ^ is a functional of a, and the action of matter (an integral over rf) is 
a functional of a. Therefore, the action cannot be separated into a sum of a gravitational part and a matter part 
as Hartle and Hawking did. The failure of Hartle and Hawking's path integral calculation is also manifested in the 
fact that de Sitter space is not a solution of the Friedmann equation if the "zero-point-energy" of the conformally 
coupled scalar field is considered, whereas the semiclassical approximation implies that the principal contribution to 
the path integral of the wave function comes from the configurations which solve Einstein's equations. One may hope 
to overcome this difficulty by introducing a scalar field with a flat potential V{4>) (as in the inflation case). But 
this does not apply to the quantum cosmology case since as a — ► the universe always becomes radiation-dominated 
unless the energy density of radiation is exactly zero (but the uncertainty principle does not allow this case to occur) . 

V. CTCS AND THE CHRONOLOGY PROTECTION CONJECTURE 

From the arguments in the last section, we find that the Universe does not seem to be created from nothing. On 
the other hand, if the Universe is created from something, that something could have been itself. Thus it is possible 
that the Universe is its own mother. In such a case, if we trace the history of the Universe backward, inevitably 
we will enter a region of CTCs. Therefore CTCs may play an important role in the creation of the Universe. It is 
interesting to note that Hawking and Penrose's singularity theorems do not apply if the Universe has had CTCs. 
And, it has been shown that, if a compact Lorentzian spacetime undergoes topology changes, there must be CTCs in 
this spacetime [ 140| , 19| , p0| . [Basically there are two type of spacetimes with CTCs: for the first type, there are CTCs 



everywhere (Godel space belongs to this type); for the second type, the CTCs are confined within some regions and 
there exists at least one region where there are no closed causal (timelike or null) curves, and the regions with CTCs 
are separated from the regions without closed causal curves by Cauchy horizons (Misner space belongs to this type). 
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In this paper, with the word "spacetimes with CTCs" we always refer to the second type unless otherwise specified.] 

While in classical general relativity there exist many solutions with CTCs, some calculations of vacuum polarization 
of quantum fields in spacetimes with CTCs indicated that the energy-momentum tensor (in this paper when we deal 
with quantum fields, with the word "the energy- momentum tensor" we always refer to "the renormalized energy- 
momentum tensor" because "the unrenormalized energy-momentum tensor" has no physical meaning) diverges as one 
approaches the Cauchy horizon separating the region with CTCs from the region without closed causal curves. This 
means that spacetimes with CTCs may be unstable against vacuum polarization since when the energy-momentum 
tensor is fed back to the semiclassical Einstein's equations (i.e. Einstein's equations with quantum corrections to the 
energy-momentum tensor of matter fields) the back-reaction may distort the spacetime geometry so strongly that a 
singularity may form and CTCs may be destroyed. Based on some of these calculations, Hawking jl9],^0| has proposed 
the chronology protection conjecture which states that the laws of physics do not allow the appearance of CTCs. (It 
should be mentioned that the chronology protection conjecture does not provide any restriction on spacetimes with 
CTCs but no Cauchy horizons since there is no any indication that this type of spacetime is unstable against vacuum 
polarization. In the next section we will show a simple example of a spacetime with CTCs but no Cauchy horizons, 
where the energy- momentum tensor is finite everywhere.) 

But, on the other hand, Li, Xu, and Liu p2| have pointed out that even if the energy-momentum tensor of vacuum 
polarization diverges at the Cauchy horizon, it does not mean that CTCs must be prevented by physical laws because: 
(1) Einstein's equations are local equations and the energy-momentum tensor may diverge only at the Cauchy horizon 
(or at the polarized hypersurfaces) and be well-behaved elsewhere within the region with CTCs; (2) the divergence 
of the energy-momentum tensor at the Cauchy horizon does not mean that the Cauchy horizon must be destroyed 
by the back-reaction of vacuum polarization, but instead means that near the Cauchy horizon the usual quantum 
field theory on a prescribed classical spacetime background cannot be used and the quantum effect of gravity must 
be considered. (This is like the case that Hawking and Penrose's singularity theorems do not mean that the Big 
Bang cosmology is wrong but mean that near the Big Bang singularity quantum gravity effects become important 
|n|.) When Hawking proposed his chronology protection conjecture, Hawking |20| and Kim and Thorne |ll| had 
a controversy over whether quantum gravity can save CTCs. Kim and Thorne claimed that quantum gravitational 
effects would cut the divergence off when an observer's proper time from crossing the Cauchy horizon was the Planck 
time, and this would only give such a small perturbation on the metric that the Cauchy horizon could not be destroyed. 
But, Hawking pQ ] noted that one would expect the quantum gravitational cut-off to occur when the invariant distance 
from the Cauchy horizon was of order the Planck length, and this would give a very strong perturbation on the metric 
so that the Cauchy horizon would be destroyed. Since there does not exist a self-consistent quantum theory of gravity 
at present, we cannot judge who (Hawking or Kim and Thorne) is right. But in any case, these arguments imply 
that in the case of a spacetime with CTCs where the energy-momentum tensor of vacuum polarization diverges at the 
Cauchy horizon, quantum gravity effects should become important near the Cauchy horizon. Li, Xu, and Liu have 
argued that if the effects of quantum gravity are considered, in a spacetime with CTCs the region with CTCs and the 
region without closed causal curves may be separated by a quantum barrier (e.g. a region where components of the 
metric have complex values) instead of a Cauchy horizon generated by closed null geodesies. By quantum processes, 
a time traveler may tunnel from the region without closed causal curves to the region with CTCs (or vice versa) , and 
the spacetime itself can also tunnel from one side to the other side of the quantum barrier p2j . In classical general 
relativity, a region with CTCs and a region without closed causal curves must be separated by a Cauchy horizon 
(compactly generated or non-compactly generated) which usually contains closed null geodesies if it is compactly 
generated p0| |. But if quantum gravity effects are considered (e.g. in quantum cosmology), they can be separated by 
a complex geometric region (as a quantum barrier) instead of a Cauchy horizon |2^ 1 . (In the path integral approach to 
quantum cosmology, complex geometries are required in order to make the path integral convergent and to ov ercome 
the difficulty that in general situations a Euclidean space cannot be directly joined to a Lorentzian space [141]). And, 
using a simple example of a space with a region with CTCs separated from a region without closed causal curves by 
a complex geometric region, Li, Xu, and Liu [^2| have shown that in such a space the energy-momentum tensor of 
vacuum polarization is finite everywhere and the chronology protection conjecture has been challenged. 

Without appeal to quantum gravity, counter-examples to the chronology protection conjecture also exist. By intro- 
ducing a spherical reflecting boundary between two mouths of a wormhole, Li |23f| has shown that with some boundary 
conditions for geodesies (e.g. the reflection boundary condition) closed null geodesies [usually the "archcriminal" for 
the divergence of the energy-momentum tensor as the Cauchy horizon is approached (see e.g. []l6|| )1 may be removed 
from the Cauchy horizon separating the region with CTCs and the region without closed causal curves. In such a 
case the spacetime contains neither closed null geodesies nor closed timelike geodesies, though it contains both closed 
timelike non-geodesic curves and closed null non-geodesic curves. Li p3| has shown that in this spacetime the energy- 
momentum tensor is finite everywhere. Following Li p3| , Low p4| has given another example of spacetime with CTCs 
but without closed causal geodesies. 

Recently, with a very general argument, Li |2q] has shown that the appearance of an absorber in a spacetime with 
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CTCs may make the spacetime stable against vacuum polarization. Li ^™ has given some examples to show that 
there exist many collision processes in high energy physics for which the total cross-sections increase (or tend to a 
constant) as the frequency of the incident waves increases. Based on these examples, Li Q has argued that material 
will become opaque for waves (particles) with extremely high frequency or energy, since in such cases the absorption 
caused by various types of scattering processes becomes very important. Based on calculation of the renormalized 
energy-momentum tensor and the fluctuation in the metric, Li |2q | has argued that if an absorbing material with 
appropriate density is introduced, vacuum polarization may be smoothed out near the Cauchy horizon so that the 
metric perturbation caused by vacuum fluctuations will be very small and a spacetime with CTCs can be stable 
against vacuum polarization. 

Boulware (2l]j and Tanaka and Hiscock [^j| have found that for sufficiently massive fields in Gott space | . 142 1 and 
Grant space p8| respectively, the energy-momentum tensor remains regular on the Cauchy horizon. Krasnikov |27| has 
found some two-dimensional spacetimes with CTCs for which the energy-momentum tensor of vacuum polarization 
is bounded on the Cauchy horizon. Sushkov |2^] has found that for an automorphic complex scalar field in Misner 
space there is a vacuum state for which the energy-momentum tensor is zero everywhere. More recently, Cassidy |3l]] 
and Li and Gott [Q have independently found that for the real conformally coupled scalar field in Misner space there 
exists a quantum state for which the energy-momentum tensor is zero everywhere. Li and Gott pof have found that 
this quantum state is the "adapted" Rindler vacuum (i.e. the usual Rindler vacuum with multiple images) and it is 
a self-consistent vacuum state because it solves the semiclassical Einstein's equations exactly. Li and Gott ||(]] have 
also found that for this "adapted" Rindler vacuum in Misner space, an inertial particle detector perceives nothing. In 
this paper, we find that for a multiply connected de Sitter space there also exists a self-consistent vacuum state for a 
conformally coupled scalar field (see section IX). 

Thorne |I has noted that, even if Hawking's argument that a quantum gravitational cut-off would occur when the 
geometric invariant distance from the Cauchy horizon is of order the Planck length is correct, by using two wormholes 
the metric fluctuations near the Cauchy horizon can be made arbitrarily small so a spacetime with CTCs created 
from two wormholes can be stable against vacuum polarization. Recently Visser [^9| has generalized this result to the 
Roman-ring case. 

The above arguments indicate that the back-reaction of vacuum polarization may not destroy the Cauchy horizon 
in spacetimes with CTCs, and thus such spac etim es can be stable against vacuum polarization. 

In a recent paper, Cassidy and Hawking 1 143 ] have adm itted that "back-reaction does not enforce chronology 
protection". On the other hand, Cassidy and Hawking [ 143 1 have argued that the "number of states" may enforce 
the chronology protection conjecture since "this quantity will always tend to zero as one tries to introduce CTCs" . 
Their arguments are based on the fact that for the particular spacetime with CTCs they constructed [which is the 
product of a multiply connected (via a boost) three-dimensional de Sitter spac e and S 1 ] the entropy of a massless 
scalar field diverges to minus infinity when the spacetime develops CTCs [143]. However, whether this conclusion 
holds for general spacetimes with CTCs remains an open question and further research is required. And, from 
ordinary statistical thermodynamics we know that entropy is always positive, so the physical meaning of a negative 
entropy is unclear. The number of states in phase space is given by TV = ApAq/(2irh) B where Aq — Aqi Aq2-..Aq s , 
Ap = ApiAp2---Ap s , qi (i — 1, 2, ...s) is a canonical coordinate, pi is a canonical momentum, and s is the number of 
degrees of freedom. The uncertainty principle demands that ApiAqi > 2nh and thus we should always have TV > 1. 
Thus the "fact" that the number of states tends to zero as one tries to develop CTCs (i.e. as one approaches the 
Cauchy horizon) may simply imply that near the Cauchy horizon quantum effects of gravity cannot be neglected, 
which is consistent with Li, Xu, and Liu's argument (^2|. The entropy is defined by fee In TV where TV is the number of 
states and &b is the Boltzmann constant. When TV is small, quantization of the entropy becomes important (remember 
that the number of states TV is always an integer). The entropy cannot continuously tend to negative infin ity; it should 
jump from A:eln3 to /cBhi2, jump from fcsln2 to zero (but in Cassidy and Hawking's arguments [143] we have not 
seen such a jump), then the uncertainty principle demands that the entropy should stand on the zero value as one 
approaches the Cauchy horizon. On the other hand, ordinary continuous thermodynamics holds only for the case 
with TV 3> 1. Thus, as one approaches the Cauchy horizon the thermodynamic limit has already been violated and 
ordin ary thermodynamics should be revised near the Cauchy horizon. In other words, Cassidy and Hawking's results 
[143 1 cannot be extended to the Cauchy horizon. Based on the fact that the effec tive action density diverges at 
the polarized hypersurfaces of spacetimes with CTCs |3lj] , Cassidy and Hawking |143| ] have argued that the effective 
action "would provide new insight into issues of chronology protection" . But we should note that the effective action 
is only a tool for computing some physical quantities (such as the energy-momentum tensor) and the effective action 
itself has not much physical meaning. The divergence of the effective action may imply that the effective action is not 
a good tool as the polarized hypersurfaces are approached. Our argument is supported by the fact that there exist 
many examples for which the energy-momentum tensor is finite everywhere, as mentioned above. 

Recently, Kay, Radzikowski, and Wald [144] have proved two theorems which demonstrate that some fundamental 
quantities such as Hadamard functions and energy-momentum tensors must be ill-defined on a compactly generated 
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Cauchy horizon in a spacetime with CTCs, as one extends the usual quantum field theory in a global hyperbolic 
spacetime to an acausal spacetime with a compactly generated Cauchy horizon. Basically spe aking , their theorems 



imply that the usual quantum field theory cannot be directly extended to a spacetime with CTCs [ 144 1 . Their theorems 



tell us that serious difficulties arise when attempting to define quantum field theory on a spacetime with a compactly 



generated Cauchy horizon [144]. The ordinary quantum field theory must be significantly changed or some new 



approach must be introduced when one tries to do quantum field theory on a spacetime with CTCs. A candidate 



procedure for overcoming this difficulty is the Euclidean quantization proposed by Hawking [ 145 : 146 1 . Quantum field 
theory is well-defined in a Euclidean space because there are no CTCs in a Euclidean space |147|| . In fact, even in 
simply connected Minkowski spacetime, quantum field theory is not well-defined since the path integral does not 
converge. To overcome this difficulty, the technique of Wick-rotation (which is essentially equivalent to Euclidean 



quantization) is used. Kay, Radzikowski, and Wald [144] have also argued that their results may be interpreted as 



indicating that in order to create CTCs it would be necessary to enter a regime where quantum effects of gravity will 



be dominant (see also the discu ssions of Visser [148 149]); this is also consistent with Li, Xu, and Liu's arguments 



144 1 also apply to Misner 



, respectively) where the 



22|. Cramer and Kay [150 151] have shown that Kay, Radzikowski, and Wald's theorems 
space (for Sushkov's automorphic field case |^8) and Krasnikov's two-dimensional case [^ 
Cauchy horizon is not compactly generated, in the sense that the energy-momentum tensor must be ill-defined on the 
Cauchy horizon itself. But we note that this only happens in a set of measure zero which does not make much sense 
in physics for if the renormalized energy-momentum tensor is zero everywhere except on a set of measure zero where 
it is formally ill-defined, then continuity would seem to require setting it to zero there also f3(i|| . 

Perhaps a conclusion on the chronology protection conjecture can only be reached after we have a quantum theory 
of gravity. However, we can conclude that the back-reaction of vacuum polarization does not enforce the chronology 
protection conjecture, a point Hawking himself also admits ]143| . (Originally the back-reaction of vacuum polarization 
was supposed to be the strongest candidate for chronology protection Jl9] , po|] .) 

VI. MULTIPLY CONNECTED MINKOWSKI SPACETIMES WITH CTCS 

A simple spacetime with CTCs is obtained from Minkowski spacetime by identifying points that are related by time 
translation. Minkowski spacetime is (R^,7] a i,). In Cartesian coordinates (t,x,y,z) the Lorentzian metric r\ a b is given 
by 

ds 2 = -dt 2 + dx 2 +dy 2 + dz 2 . (10) 

Now we identify points (t,x,y,z) with points (f + nto, x,y, z) where to is a positive constant and n is any integer. 
Then we obtain a spacetime with topology S 1 x R 3 and the Lorentzian metric. Such a spacetime is closed in the time 
direction and has no Cauchy horizon. All events in this spacetime are threaded by CTCs. (This is the only acausal 
spacetime without a Cauchy horizon considered in this paper.) Minkowski spacetime (i? 4 ,?7 a fc) is the covering space 
of this spacetime. 

Usually there is no well-defined quantum field theory in a spacetime with CTCs. (Kay-Radzikowski- Wald's theorems 



[144 1 enforce this claim, though they do not apply directly to an acausal spacetime without a Cauchy horizon.) 
However, in the case where a covering space exists, we can do it in the covering space with the method of images. In 
fact in most cases where the energy-momentum tensor in spacetimes with CTCs has been calculated, this method has 
been used (for the theoretical basis for the method of images see Ref. and references cited therein) . The method of 
images is sufficient for our purposes in this paper (computing the energy-momentum tensor and the response function 
of particle detectors). Thus in this paper we use this method to deal with quantum field theory in spacetimes with 
CTCs. 

For any point (£, x, y, z) in (S 1 x i? 3 , ?7ab), there are an infinite number of images of points (t + nto, y, z) in the 
covering space (i? 4 , 77^5) . For the Minkowski vacuum |0m) of a conformally coupled scalar field (by "conformally 
coupled" we mean that the mass of the scalar field is zero and the coupling between the scalar field <fi and the 
gravitational field is given by \R(j) 2 where R is the Ricci scalar curvature) in the Minkowski spacetime, the Hadamard 
function is 

r<( 1 ) (y x'\ = - - f 1 1*1 

mK ' ' 2ir 2 ~(t-t>) 2 + (x-x') 2 + (y-y') 2 + (z- z') 21 1 ' 

here X = (t,x,y,z) and X' = (t', x', y', z'). With the method of images, the Hadamard function of the "adapted" 
Minkowski vacuum (which is the Minkowski vacuum with multiple images) in the spacetime (S 1 x R 3 ,r] a b) is given 
by the summation of the Hadamard function in (Fyl) for all images 
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1 oo 



2tt 2 ^-^ - t' + nt Q ) 2 + (x - x') 2 + (y - y') 2 + (z - z') 2 ' 



(12) 



The regularized Hadamard function is usually taken to be 
g£Ux,x>) = G^(X,X>) - G$(X t X') 



—y 



2tt 2 ^ -{t -t' + nt ) 2 + {x- x 1 ) 2 + (y - y') 2 + (z - z') 2 ' 



The renormalized energy-momentum tensor is given by [152|]153[] 



(T ab ) 

Inserting Eq. ( |l3| ) into Eq. (|l4|) we get 



— lim 

2 x'^x 



;V a V fc 



1. 



1 



90^ 



;V a V b - -?? ah V c V 
6 



-3 
10 
10 
1 



(13) 



(14) 



(15) 



We find that this energy-momentum tensor is constant and finite everywhere and has the form of radiation. Thus 
CTCs do not mean that the energy-momentum tensor must diverge. 

Now let us consider a particle detector [152,154] moving in this spacetime. The particle detector is coupled to 
the field by the interaction Lagrangian crn[T)^[X (t)] , where c is a small coupling constant, m is the detector's 
mon opole moment, r is the proper time of the detector's worldline, and X(t) is the trajectory of the particle detector 
[ 153 1 . Suppose initially the detector is in its ground state with energy Eq and the field (f> is in some quantum state |). 
Then the transition probability for the dete ctor to all possible excited states with energy E > Eq and the field </> to 
all possible quantum states is given by |153 || 

(16) 



P = c 2 £ K J E|m(0)|E )| 2 ^(A^) ) 

E>E 



where AE = E — Eq > and T(AE) is the response function 

dr / dT'e- lAE(T - T ">G+(X(T),X(T')), 

-oo J — oo 



(17) 



which is independent of the details of the particle detector and is determined by the positive frequency Wight- 
man function G+(X,X') = (\<p(X)(f>(X')\) (while the Hadamard function is defined by G 1 (X, X') = (\4>(X)<f)(X') + 
4>(X')(f>(X)\)). The response function represents the bath of particles that the detector effectively experiences [ 153 1 . 
The remaining factor in Eq. ( |l6| ) represents the selectivity of the detector to the field and depends on the internal 
structure of the detector [153]. The Wightman function for the Minkowski vacuum is 



@m(X, X') 



4tt 2 -(t - t' - ie) 2 + (x - x') 2 + (y - y') 2 + (z - z') 2 ' 



(18) 



where e is an infinitesimal positive real number which is introduced to indicate that G + is the boundary value of a 
function which is analytic in the lower-half of the complex At = t — t' plane. For the adapted Minkowski vacuum in 
our spacetime (S 1 x R 3 ,r] a b), the Wightman function is 



G + (X,X') 



1 

47T 2 



E 



(19) 



n— — oo 



Assume that the detector moves along the geodesic x = (3t ((3 < I), y = z 
C = 1 / \J 1 — [i 2 . On the geodesic, the Wightman function is reduced to 



-(t - t' + ntQ - it) 2 + (x - x') 2 + (y- y') 2 + (z - z') 2 ' 

0, then the proper time is r = t/ Q with 
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(t-<' + 7rf o -«e) 2 + /3 2 (t-<0 2 

1 00 1 

£ ( T - r ' + n i /C - «e/C) 2 - P 2 (r - r') 2 ' (20) 



n— — 00 



Inserting Eq. ( pp| ) into Eq. Jl7|), we obtain 

1 °° p 00 /» 00 -1 
F(AE) = - — - dT dAre- iAEA T - — — — -. (21) 

4?r C 2 „ J-oc J -00 ( Ar + « VC - WO 2 - /? 2 (At) 2 



n— — 00 



where Ar = r — r' and T = (r + r')/2. The integration over At is taken along a contour closed in the lower-half plane 
of complex At. Inspecting the poles of the integrand, we find that all poles are in the upper-half plane of complex 
At (remember that [3 < 1). Therefore according to the residue theorem we have 

T{AE) = 0. (22) 

Such a particle detector perceives no particles, though the renormalizcd energy-momentum tensor of the field has the 
form of radiation. 

Another simple space with CTCs constructed from Minkowski space is Misner space ||. In Cartesian coordinates 
(t, x, y, z) in Minkowski spacetime, a boost transformation in the (t, x) plane (we can always adjust the coordinates 
so that the boost is in this plane) takes point (f , x,y, z) to point (t cosh b + x sinh b, x cosh b + t sinh b, y, z) where b is 
the boost parameter. In Rindler coordinates (f],^,y, z), defined by 



t = £ sinh 77, 
x = £ cosh r], 

y = y, 

z — z, 



(23) 



the Minkowski metric can then be written in the Rindler form 

ds 2 = -fdrf + d£ 2 + dy 2 + dz 2 . (24) 

The Rindler coordinates (r),£,y,z) only cover the right quadrant of Minkowski space (i.e. the region R defined by 
x > \t\). By a reflection (t,x,y,z) — > (—t,—x,y,z) [or (rj,£,y, z) — > (17, — £, y, z)] , the Rindler coordinates and the 
Rindler metric can be extended to the left quadrant (L, defined by x < — \t\). By the transformation 

7T 

»7 £ -> ±ifj, y^y, z -> z, (25) 

the Rindler coordinates can be extended to the future quadrant (F, defined by t > |x|) and the past quadrant (P, 
defined by t < — \x\). In region L the Rindler metric has the same form as the metric in region R, which is given by 
Eq. (p3). But in F and P the Rindler metric is extended to be 

ds 2 = -dfj 2 + tfdti 2 + dy 2 + dz 2 . (26) 

Misner space is obtained by identifying (t, x,y, z) with (t cosh nb + x sinh nb, x cosh nb + t sinh nb, y, z) . Under such an 
identification, point (r),^,y,z) in R (or L) is identified with points (77 + nb, £, y, z) in R (or L), point (77, £, y, z) in F 
(or P) is identified with points (fj, £ + nb, y, z) in F (or P). Clearly there are CTCs in R and L but there are no closed 
causal curves in F and P, and these regions are separated by the Cauchy horizons x — ±t, generated by closed null 
geodesies. 

Misner space is not a manifold at the intersection of x — t and x = —t. However, as Hawking and Ellis have 
pointed out, if we consider the bundle of linear frames over Minkowski space, the corresponding induced bundle of 
linear frames over Misner space is a Hausdorff manifold and therefore well-behaved everywhere. 

The energy-momentum tensor of a conformally coupled scalar field in Misner space has been studied in Jl4] , p0[ . 
Hiscock and Konkowski fl4|| have calculated the energy-momentum tensor of the adapted Minkowski vacuum. In 
Rindler coordinates their results can be written as 



^4 

(T u u )m,tci 



12tt 2 £ 4 



/ -3 
10 
10 
1 



(27) 
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where the constant A is 



a _ 2 + coshnfr 
A -Z^ ( C o S hn&-l) 2 - 



n=l 



Eq. (P7| ) holds only in region R [because Rindler coordinates defined by Eq. (g3|) only cover R] , but it can be analytically 
extended to other regions by writing (T^M.rcn in Cartesian coordinates or by the transformations mentioned above. 
Obviously for any finite b, (X^^M.rcn diverges as one approaches the Cauchy horizon (£ — > 0). This divergence is 
coordinate independent since (T^M.renC^u^M.ren a l so diverges as £ — > 0. This indicates that though the Minkowski 
vacuum is a good and self-consistent vacuum for simply connected Minkowski space, the adapted Minkowski vacuum 
is not self-consistent for Misner space (i.e. it does not solve Einstein's equations given the Misner space geometry). 
This result has led Hawking Jl9],^0| 'o conjecture that the laws of physics do not allow the appearance of CTCs (i.e., 
his chronology protection conjecture). 

Li and Gott |30[ have studied the adapted Rindler vacuum in Misner space. The Hadamard function for the Rindler 



vacuum is [155 



where X = (rj, £, y, z), X' — (77', y', z'), and 7 is defined by 

cos]i7 _ e + e* + (y- yY + U-t')\ (30) 

The Hadamard function for the adapted Rindler vacuum in Misner space is 

1 00 7 

G(1) X>S) = 2^2 sinh7[-(7 ? -r)'+ nbf + 7 2 ] ' (31) 



Though G^ and G (1) given by Eq. (^9|) and Eq. (|3l]) are defined only in region R, they can be analytically extended 
to regions L, F, and P in Minkowski and Misner space. The regularized Hadamard function for the adapted Rindler 
vacuum is Grll(X, X') = G^ {X, X') — G^ {X, X'), where G^ is the Hadamard function for the Minkowski vacuum 
given by Eq. ([ll]). Inserting this together with Eq. |3l]) and Eq. ([ll]) into Eq. (H), we obtain the energy- momentum 
tensor for a conformally coupled scalar field in the adapted Rindler vacuum |3Q] 



(T/> R , rcn - M407r2e4 

















(t) - 1 
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(32) 



which is expressed in Rindler coordinates and thus holds only in region R but can be analytically extended to other 
regions with the method mentioned above for the case of the adapted Minkowski vacuum. We ]3(| have found that 
unless b = 2ir, (T^)^ Icn blows up as one approaches the Cauchy horizon (£ — ► 0) (as also does(T M1/ )ri iren (T Atl/ )R, rGn ). 
But, if 6 = 2ir, we have 

(T/) R , rea = 0, (33) 

which is regular as one approaches the Cauchy horizon and can be regularly extended to the whole Misner space, 
where it is also zero. In such a case, the vacuum Einstein's equations without cosmological constant are automatically 
satisfied. Thus this is an example of a spacetime with CTCs at the semiclassical quantum gravity level. We [^o| have 
called this vacuum the self-consistent vacuum for Misner space, and b = 2n is the self-consistent condition. (Cassidy 
ifjlf has also independently proven that for a conformally coupled scalar field in Misner space there should exist a 
quantum state for which the energy-momentum tensor is zero everywhere. But he has not shown what quantum state 
it should be. We ]3(| have shown that it is the adapted Rindler vacuum.) 

Another way to deal with quantum fields in spacetimes with CTCs is to do the quan tum field theory in the 
Euclidean section and then analytically extend the results to the Lorentzian section ]147[ . For Misner space the 
Euclidean section is obtained by taking 77 and b to be —if) and — ib. The resultant space is the Euclidean space 
with metric ds 2 = £ 2 dfj 2 + d£ 2 + dy 2 + dz 2 and {fj,£,y,z) and (fj + nb,£,y,z) are identified where (fj,£,y,z) are 
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cylindrical polar coordinates with fj the angular polar coordinate and £ the radial polar coordinate. The geometry 
at the hypersurface £ = is conical singular unless b = 2ix. When extending that case to the Lorentzian section, we 
get b = 2tt which is just the self-consistent condition. This may be the geometrical explanation of the self-consistent 
condition. By doing quantum field theory in the Euclidean space, then analytically extending the results to the 
Lorentzian section, we obtain the renormalized energy-momentum tensor in R (or L) region of the Misner space. 
Then we can extend the renormalized energy- momentum tensor in R (or L) to regions F (or P). The results are the 
same as that obtained with the method of images. 

Let us consider a particle detector moving in Misner space with the adapted Rindler vacuum. Suppose the detector 
moves along a geodesic with x = a, y = (3t, and z = (a and f3 are constants and a is positive), which goes through 
the P, R, and F regions. The proper time of the detector is r = £/£ with Q = 1/yl — /3 2 . On this geodesic, the 
Hadamard function in (31) is reduced to 



1 



7 



oo 



where 7 is given by 



2vr 2 



smh-y^J (a 2 -t 2 )(a 2 -t' 2 ) n=-oo 



1 



-(77 — rf + nb) 2 + 7 2 



(34) 



cosh 7 



2a 2 -t 2 + f3 2 {t-t') 



A2 



2J(a 2 -t 2 )(a 2 -t' z ) 



and r) — rj is given by 



sinh(?7 — rj ) 



a(t - t') 
yj '(a 2 - t 2 ){a 2 - t' 2 ) 



(35) 



(36) 



Though this Hadamard function is originally defined only in R, it can be analytically extended to F, P, and L. The 
Wightman function is equal to 1/2 of the Hadamard function with t replaced by t — ie/2 and t' replaced by t' + ie/2, 
where e is an infinitesimal positive real number. Then the response function is |30(] 



-1 x roc roc 

^-'—oo J— oc 



7 ^e 



— iEAr 



sinl^+^/K - C 2 (T + *f - ||) 2 ][a 2 - C 2 (T - $ + %) 2 ] {-[fa-Tj'H 



ib] 2 + 7 ^ 



(37) 



where T = (r + r')/2, At = r — r'; 7 + and (77 — rj') + are given by ( p5[) and (36) with i replaced by i — ie/2 and t' 



replaced by t' + ie/2. The integral over At can be worked out by the residue theorem where we choose the integration 
contour to close in the lower-half complex-Ar plane. The result is zero since there are no poles in the lower-half plane. 
Therefore such a detector cannot be excited and so it detects nothing |50|. We p(J have also calculated the response 
functions for detectors on worldlines with constant £, y, and z and worldlines with constant £, y, and z — both are 
zero. 



VII. VACUUM POLARIZATION IN VILENKIN'S TUNNELING UNIVERSE 



In order to compare our model for the creation of the universe with Vilenkin's tunneling universe, in this section we 
calculate the vacuum fluctuation of a conformally coupled scalar field in Vilenkin's tunneling universe. The geometry 
of Vilenkin's tunneling universe has been described in section IV. Such a universe is described by a Lorentzian-de Sitter 
space joined to a Euclidean de Sitter space |ll[] . The Lorentzian section has the topology R 1 x S 3 and the metric 

ds 2 = -dr 2 + r 2 cosh 2 y[d X 2 + sin 2 X (d0 2 + sin 2 9d(j) 2 )]. (38) 

The Euclidean section has the topology S 4 and the metric 

ds 2 = dr 2 + r 2 cos 2 — [d X 2 + sin 2 X {d0 2 + sin 2 8d<p 2 )}. (39) 
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The Lorentzian section and the Euclidean section are joined at the boundary S defined by r = 0. S is a three-sphere 
with the minimum radius in de Sitter space and th e maxim um radius in the Euclidean four-sphere. The boundary 
condition for a conformally coupled scalar field <f> is [156,157] 



dr 



0. 



(40) 



which is a kind of Neumann boundary condition and indicates that the boundary E is like a kind of reflecting boundary. 
The Green functions (including both the Hadamard function and the Wightman function) should also satisfy this 
boundary condition 



0G(t 1 x,M;t , ,x / ,0',^) 



0. 



(41) 



The vacuum state of a conformally coupled scalar field in de Sitter space is usually taken to be that obtained from 
the Minkowski vacuum by the conformal transformation according t o wh ich de Sitter space is conformally flat. (The 
quantum state so obtained is usually called the conformal vacuum [153].) Such a vacuum is de Sitter invariant and 
we call it the conforma l Mi nkowski vacuum. The Hadamard function for this de Sitter vacuum (i.e. the conformal 
Minkowski vacuum) is [158] 



1 



1 



where X = (r, Xi 



A-K 2 rl 1-Z(X,X')' 
X' = (r', x', <?', ft), and Z(X, X') is defined by 



(42) 



T T T T 

Z(X, X') = — sinh — sinh h cosh — cosh — {cos x cos x' 

ro r r r 

+ sin x sin x' [cos 9 cos 9' + sin sin cos(0 — </>')]}■ 
In Vilcnkin's tunneling universe, the Hadamard function satisfying the boundary condition ([il"| ) is given by 

G« (X, X') = G<& (X, X') + G« (X- , X') 
1 



(43) 



47T 2 ? 



1-Z(X,X') 1-Z(X-,X') 



(44) 



where X = (— r, x, 6, 4>) is the image of X — (r, x, 9, <f>) with respect to the reflecting boundary E. 

There are various schemes for obtaining the renormalized energy-momentum tensor for de Sitter space (e.g. 

[158,159]). They all are equivalent to subtracting from the Hadamard function a reference term GjH to o btain a 
regularized Hadamard function and then calculating the renormalized energy-momentum tensor by [152 153 



(Tab 



±Km x V ab ,(X,X')G£l(X,X') 



For the conformally coupled scalar field, the differential operator T) ao i is 



v aV = -v a v fc - 
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1 



1 



77 Rab 
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2 R 9ab 



(45) 



(46) 



whereoafc' is the geodesic parallel displacement bivector [16C|. [It is easy to show that if R a b = Eq. ( [45[ ) and 
Eq. ( [iq ) are reduced to Eq. ( |l4| ) .] The regularized Hadamard function for the adapted conformal Minkowski vacuum 
in Vilenkin's tunneling universe is 



Ggl(X,X') = GU(X,X') G { V{X,X') = G« (X,X') - Gg> 



Gcm(X ,X'). 



(47) 



(In this paper the exact form of G^] is not important for us.) Substituting Eqs. ([I2j|44|) and Eq. ([47j) into Eq. (|4 
we find that limx'->x T> ^iGq^X - , X') = 0, which shows that the boundary condition ([icj ) does not produce any 
renormalized energy-momentum tensor; but the action of T> a y on 

4m(^, X')-?^ should give the energy-momentum 
tensor for the conformal Minkowski vacuum in an eternal de Sitter space [ 158) , 159( ] 
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- lim Vav 

2 x'^x 



^CMl^' A ) ~ °rcf 



(1) 



9607r 2 r4 



9ab- 



(48) 



Therefore, the energy-momentum tensor of a conformally coupled scalar field in the adapted Minkowski vacuum in 
Vilenkin's tunneling universe is 



(Ta 



1 



6/rcn 



9607r 2 r-4 



9ab, 



(49) 



which is the same as that for an eternal de Sitter space. 

Now consider a particle detector moving along a geodesic with Xi 9, 4> = constants. The response function is given 
by Eq. (17) but with the integration over r and r' ranging from to oo. The Wightman function is obtained from 
the corresponding Hadamard function by the relation 



G+(t, x, 0, </>; r', X ', 8', 0') = (t - ie/2, X , 6, 0; r' + ie/2, X ', 9', <j>') , 

where e is an infinitesimal positive real number. Along the worldline of the detector, we have 



Z(r, t') = — sinh — sinh — + cosh — cosh — 



cosh ■ 



ro 



'o 



'"o 



'0 



T + T' 



T . T T T 

sinh — sinh — + cosh — cosh — = cosh 



and 



(50) 

(51) 
(52) 



G + {X,X') 
Then the response function is 



1 - cosh l^-L-^L 



1 - cosh 



1 /• OC f- oc 



-iAEr Ar 



1 



1 - cosh( At - ie) 1 - cosh 2T 



(53) 



(54) 



where At = (t — r')/ro and T = (r + r')/2ro. It is easy to calculate the contour integral over At. We find that the 
integration of the second t erm is ze ro and therefore, the result is the same as that for an inertial particle detector in 
an eternal de Sitter space |l00|JT5^ 1. Thus we have 



dT 



AE 



2ir e 27Tr ° AE - 1 : 



(55) 



which is just the response function for a detector in a thermal radiation with the Gibbons-Hawking temperature [IOC] 

1 



G-H 



2nr 



(56) 



[The factor ro over 2-7T in Eq. ( |55| ) is due to the fact that by definition T — (r + r')/2ro is dimensionless.] Therefore 
such a detector perceives a thermal bath of radiation with the temperature Tq_h- 

Though the boundary between the Lorentzian section and the Euclidean section behaves as a reflecting boundary, 
a particle detector cannot distinguish Vilenkin's tunneling universe from an eternal de Sitter space, and they have 
the same energy-momentum tensor for the conformally coupled scalar field. 



VIII. A TIME-NONORIENTABLE DE SITTER SPACE 



A time-nonorientable de Sitter space can be constructed from de Sitter space by identifying antipodal points |161,io[. 
Under s uch an id entification, point X = (r, x, d, (f>) is identified with —X = (— t, tt — x, tt — 9, -k + </>). Friedman and 
Higuchi [162,163] have described this space as a "Lorentzian universe from nothing" (without any Euclidean section), 
although one could also describe it as always existing. Friedman and Higuchi hav e studied quantum field theory in 
this space but have not calculated the renormalized energy-momentum tensor |l62j] . 
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Dc Sitter space is the covering space of this time-nonorientable model. Using the method of images, the Hadamard 
function of a conformally coupled scalar field in the time-nonorientable de Sitter space with the "adapted" conformal 
Minkowski vacuum can be constructed as 



gU(x,x') = g£1(x,x') + g£I l (-x,x') = 



(1) 



1 



AlT 2 7 



1 



1 



1-Z(X,X') 1-Z(-X,X') 



l-Z(X,X') 1 + Z(X,X') 



(57) 



The regularized Hadamard function is 

g£Ux,x') = G^(X,X') - G^(X,X') 



Gcm(- X ' X ') ~ G^JiX^X') + G^i^—X, X') 



(58) 



Inserting Eq. 
tensor is zero. T 



and Eq. ([38]) into Eq. (ff5|), we find that the contribution of G^^— X, X') to the energy- momentum 
lerefore the renormalized energy-momentum tensor is the same as that in an eternal de Sitter space, 
which is given by Eq. (|49|). 

Suppose a particle detector moves along a worldline with x,6,4> — constants. The response function is given by 
Eq. ([L7j) . The Wightman function is obtained from the Hadamard function through Eq. ( pO| ) . On the worldline of the 
particle detector, we have 



G+(r,r') 



1 



1 



f 



Inserting this into Eq. (17) we get 



8vr 2 ^ \ 1 - cosh i^L^n i + cos h 



dT r AE 



r — r'—ia 



dT 2tt e^oAB _ i 



(59) 



(60) 



which represents a thermal spectrum with a temperature equal to twice the Gibbons-Hawking temperature. Therefore 
a particle detector moving along such a geodesic in this time-nonorientable spacetime perceives thermal radiation with 
temperature T — 2Tq-h- 

For this time-nonorientable de Sitter space, the area of the event horizon is one half that of an eternal de Sitter 
space. This together with T = 2Tg_h tells us that the first thermodynamic law of even t horizons 5M C = TSA is 
preserved, where M c is the mass within the horizon, and A is the area of the horizon [ 100 1 . 



IX. A MULTIPLY CONNECTED DE SITTER SPACE WITH CTCS 

A. Construction of a Multiply Connected de Sitter Space 

De Sitter space is a solution of the vacuum Einstein's equations with a positive cosmological constant A, which 
is one of the maximally symmetric spacetimes (the others being Minkowski space and anti-de Sitter space) p^JfOf . 
De Sitter space can be represented by a timelike hyperbolic hypersurface 

W 2 +X 2 +Y 2 + Z 2 - V 2 = r%, (61) 

embedded in a five-dimensional Minkowski space (V, W, X, Y, Z) with the metric 

ds 2 = -dV 2 + dW 2 +dX 2 +dY 2 +dZ 2 , (62) 

where r = (3/ A) 1 / 2 De Sitter space has ten killing vectors — four of them are boosts, and the other six are 

rotations. The global coordinates (r, %, 6, <f>) have been described in previous sections. Static coordinates (t,r,8,<f>) 
on de Sitter space are defined by 

' U = (r2-r 2 ) 1 /2 sin h^, 
W = (r 2 -r 2 ) 1 / 2 cosh^, 
< X = r sin 6 cos 0, (63) 
Y = r sin 9 sin 0, 
Z = r cos 0, 
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where — oo <i<oo,0<r<ro,0<#<7r, and < <f> < 2n. In these coordinates the de Sitter metric is written as 



ds 2 



— ) dt 2 



I-3 dr 2 +r 2 (d9 



2 1 sin 6 2 c 



We divide de Sitter space dS into four regions 



K = {pe dS\W > \V\}, 
£ = { P EdS\W < -\V\}, 
T={p<^dS\V > \W\}, 
P = {pedS\V < -\W\}, 



(64) 



(65) 
(66) 
(67) 
(68) 



which are separated by horizons where W — ±V and X 2 + Y 2 + Z 2 = r 2 ,. (See Fig. |^). It is obvious that the 
static coordinates defined by Eq. ( |63[ ) only cover region 1Z. However, similar to the Rindler coordinates, these static 
coordinates can be extended to region T by the complex transformation 



t -> I 



t, 



(69) 



where —00 < I < 00 and t > 2tq. In region J 7 , with the coordinates (t, I, 9, 4>), the de Sitter metric can be written as 



ds 2 



- 1 



dt 1 



-l)dl 



P{d9 2 



sm f 



Transforming the coordinate t to the proper time r by 



t = tq cosh ■ 



the de Sitter metric in T is written as 

ds 2 = -dr 2 



ro 



sinh 2 — dl 2 + rl cosh 2 — (d9 2 



sm 



(70) 



(71) 



(72) 



(See Fig. EL) The coordinates (r, /, 9, 0) are related to (V, W, X, Y, Z) by 



' V — rn sinh — cosh — , 
W = rn sinh — sinh — , 
X = rn cosh — sin 9 cos <b, 
Y = ro cosh — sin 9 sin 6, 
Z = rn cosh — cos 9. 

v u rn 



(73) 



The universe with metric ( p2\j is a type of Kantowski-Sachs universe [164|. Any hypersurface of r = constant has 
topology R 1 x S 2 and has four killing vectors. Similarly, the static coordinates can also be extended to V and C. 
Another coordinate system which will be used in this paper is the steady-state coordinate system (t, x, y, z), defined 

by 



t = rn In 



_ r X 

W+V ' 
_ rpY 
~ W+V ' 
— r Z 
W+V 



W+V 
ro ' 



(74) 



These coordinates cover regions 1Z + T and the horizon at W = V > 0. With these steady-state coordinates, the 
de Sitter metric can be written in the steady-state form 



ds 2 



-dr 2 + e 2T ^(dx 2 + dy 2 + dz 2 ). 



(75) 



Introducing the conformal time 



V = -roe 



-r/ra 



W + V 



(76) 
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and spherical coordinates (p,9,<p) denned by x — psin#cos0, y = psin^sin^, and z = pcosO, the de Sitter metric 
can be written as 

2 

ds 2 = § \-dfj 2 + dp 2 + p 2 {d9 2 + sin 2 9d0 2 )} . (77) 

Tj 

The dc Sitter metric is invariant under the action of the de Sitter group. Because the boost group in de Sitter space 
is a sub-group of the de Sitter group, the de Sitter metric is also invariant under the action of the boost group. A boost 
transformation in the (V, W) plane in the embedding five-dimensional Minkowski space induces a boost transformation 
in the de Sitter space. Under such a transformation, point (V, W, X, Y, Z) is taken to (V cosh b + W sinh b, W cosh b + 
U sinh 6, X, Y, Z). In static coordinates in 1Z, point (t, r, 9, <fi) is taken to (t + j3, r, 9, <fi) where (3 — br$. In coordinates 
(t, /, 9, (j>) in T, point (t, I, 9, <j>) is taken to (t, I + (3,9, (p). Similar to Misner space, our multiply connected de Sitter 
space is constructed by identifying points (V, W, X, Y, Z) with ( V cosh nb + W sinh nb, W cosh nb+V sinh nb, X, Y, Z) 
on de Sitter space dS. In regions 1Z, points (t, r, 9, 4>) are identified with (t+n(3, r, 9, <f>); in region T , points (t, 1, 9, <f>) are 
identified with (t, l+n(3, 9, </>). We denote the multiply connected de Sitter space so obtained by dS/B, where B denotes 
the boost group. Under the identification generated by the boost transformation, clearly dS/B has CTCs in regions 
TZ and C, but has no closed causal curves in regions T and V . The boundaries at W = ±V and X 2 + Y 2 + Z 2 = r 2 , are 
the Cauchy horizons which separate the causal regions T and V from the acausal regions 1Z and C and are generated 
by closed null geodesies (Fig. ||). 

Similar to the case of Misner space, dS/B is not a manifold at the two-sphere defined by W = V = and 
X 2 +Y 2 + Z 2 = Tq. However, as in Hawking and Ellis's arguments for Misner space p0| , the quotient of the bundle 
of linear frames over de Sitter space by the boost group is a Hausdorff manifold and thus is well-behaved everywhere. 
It may not be a serious problem in physics that dS/B is not a manifold at the two-sphere mentioned above since this 
is a set of measure zero. 



B. Conformal Relation between Our Multiply Connected de Sitter Space and Misner Space 



It is well known that de Sitter space is conformally flat, 
by the conformal transformation 



The de Sitter metric is related to the Minkowski metric 



Sab = ^ 2 Vab- 



(78) 



It is easy to show this relation by writing the steady-state de Sitter metric using conformal time [see Eq. ([n])]. 
However, in this paper it is more convenient to show this conformal relation by writing the de Sitter metric in the 
static form and the Minkowski metric in the Rindler form, and using the transformation [ 165 1 



then the conformal factor fl 2 is 



fi 2 



r > = r > 



_ y/l-r 2 /rl 



1—7' cos 0/ro 5 

r sin 8 cos 0/ro 

y l — r cos 6 /ro ' 

r sin sin 0/ro 

l—r cos 0/ro ' 



Tq(1 - rcos0/r ) 2 



(79) 



(80) 



The conformal relations given by Eq. (|79|) and Eq. ( j80| ) define a conformal map between the static de Sitter space and 
the Rindler space. The horizon at r = ro in the static de Sitter space coordinates corresponds to the horizon £ = 
in Rindler space, and the worldline r = in de Sitter space corresponds to the worldlinc with £ = 1 and y — z = in 
Rindler space. This conformal relation can also be extended to region T in de Sitter space and region F in Minkowski 
space, where we have 



1—t cos B fro 7 



r 7 

t sin cos 0/ro 
1 — t cos /ro : 

t sin 6 sin 0/ro 
1 — t cos 0/ro ' 



(81) 
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and 

Q 2 = r 2 (l - tcos9/r ) 2 . (82) 



Eq. (81) and Eq. (82) give a locally conformal map in the sense that in T in de Sitter space, the map give n by Eq. (81) 
and Eq. ( |82] ) with a "+" sign only covers 9 < 9 < it, where 8 — Arccos(r /i); the map given by Eq. ( |Sl| ) and Eq. (82) 
with a "— " sign only covers < 9 < 9q. (Remember that in F in Rindler space we have ij > 0.) This conformal 
map is singular at 9 = 9q. However, since the hypersurfaces t = constant and ij = constant are homogeneous, in a 
neighborhood of any point in region F, we can always adjust coordinates (9, (f>) so that Eq. ( pi] ) and Eq. ( |82| ) hold, 
except for the points lying in region O defined by ij 2 > 1 + y 2 + z 2 (i.e. t 2 — x 2 — y 2 — z 2 > 1) in F; because as t — ► oo 
we have i) 2 / (1 + y 2 + z 2 ) — > 1. This means that there always exists a locally conformal map between T and F-0 
(defined by t 2 — x 2 — y 2 — z 2 < 1 in F), and future infinity (t — > oo) in T corresponds to the hyperbola rj 2 = 1 + y 2 + z 2 
(i.e. t 2 — x 2 — y 2 — z 2 = 1) in F. 

With the above conformal transformation, Misner space is naturally transformed to the multiply connected de Sitter 
space dS/B with 

p = br . (83) 

For a conformally coupled scalar field in a conformally flat spacetime, the Green function G(X , X') of the conformal 
vacuum is related to the corresponding Green function G(X, X') in the flat spacetime by [ 153 1 

G{X,X') = QT l {X)G{X,X')Sr 1 {X'), (84) 

the renormalized energy-momentum tensors are related by ]153| 



(T a b ) ien = Sl- i (T a b ) 



16tt 2 
where 



(85) 



^H ab = 2V a V 6 i? - 2g ab V c V c R - ^R 2 g ab + 2R,R ab , (86) 

2 1 1 

^H ab = R^R cb — -RR ab — -R c dR cd g a b + -^R 2 g a b, (87) 



and for scalar field we have a\ = and 03 = —3^ |l53fl . [The sign before 1/167T 2 is positive here because 
we are using signature (—,+,+,+)]. For de Sitter space we have R ab — Ag ab , R — 4A, and thus ^' H ab = 0, 
^H ab = ^A 2 g ab = ^g a b- Inserting them into Eq. ([35]), we have 



7= b, 



(T a ») Icn = n-*(T a ) r 



en 



9607r 2 r4 Q 



Since the renormalized energy-momentum tensor for Minkowski space in the Minkowski vacuum is zero, we have 

(T a ) rcn = 0, and thus for a conformally coupled scalar field in the conformal Minkowski vacuum in a simply connected 
de Sitter space dS 

(Tab/reii = ~ T^JT; n~T9ab, (89) 

which is just the expected result [see Eq. (^9j)]. 

If we insert the energy-momentum tensor in Eq. ( |89| ) into the semiclassical Einstein's equations 

G ah + A. gab = 87r(T ab ) ron , (90) 

and recall that for de Sitter space we have G ab = R ab — hRgab = —-^zQab, we find that the semiclassical Einstein's 
equations are satisfied if and only if 

A -^ + W = - (91) 
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If A = 0, the solutions to Eq. (91) are vq — (3607r) -1 / 2 and ro — oo J89|] . Gott [p9| has called the vacuum state in 
de Sitter space wit h rp = (3607r) / 2 the self-consistent vacuum state (it has a Gibbons-Hawking thermal temperature 
Tg-h = l/27rr ) ItLOOj . In this self-consistent case, (T Q (,) ren = — <7 a b/9607r 2 7-Q itself is the source term producing the 
de Sitter geo metry |89fl . This may give rise to inflation at the Planck scale J89| . (In a recent paper of Panagiotakopoulos 
and Tetradis [167], inflation at the Planck scale has been suggested to lead to homogeneous initial conditions for a 
second stage inflation at the GUT scale.) The second solution ro = oo corresponds to Minkowski space. These 
perhaps supply a possible reason that the effective cosmological constant is either of order unity in Planck units or 
exactly zero. That is interesting because we observe A c ff = today and a high A c g is needed for inflation. If A ^ 0, 
we find that the solutions to Eq. (J9l|) are 



lit 1 



A 



270vr 



(92) 



A de Sitter space with ro given by Eq. (^2[) automatically satisfies the semiclassical Einstein's equations (|90|). Such a 
de Sitter space and its corresponding vacuum are thus self-consistent. 



C. Renormalized Energy-Momentum Tensor in Multiply Connected de Sitter Space 



From Eq. (|8S]) we find that if we know the energy-momentum tensor of a conformally coupled scalar field in some 
vacuum state in Misner space, we can get the energy-momentum tensor in the corresponding conformal vacuum in 
the multiply connected de Sitter space. 



Two fundamental vacuums in Minkowski space are the Minkowski vacuum and the Rindler vacuum [152, 16£]. The 
energy-momentum tensor of the conformally coupled scalar field in the adapted Minkowski vacuum in Misner space 
has been worked out by Hiscock and Konkowski juj; their results are given by Eq. Inserting Eq. ( |27| ) into 

Eq. (pH), and using Eqs. (|79|-p3[), we obtain the energy-momentum tensor of a conformally coupled scalar field in 
the adapted conformal Minkowski vacuum in our multiply connected de Sitter space dS/B. In static coordinates 
(t, r, 6, cj)), it is written as 



A 



CM.rcn 



12tt 2 4 (1 



r 2 /r 2 ) 2 



-3 
10 
10 
1 



960tt 



2-4 



s v 



(93) 



where 



A = 



S (cosh 22 -l)' 



cosh 2fi 

ro 



(94) 



This result is defined in region 7Z, but it can be extended to region T through the transformation in Eq. (|69|) , and can 
also be extended to region C and V through similar transformations. Similar to Misner space, this energy-momentum 
tensor diverges at the Cauchy horizon as r — > ro for any finite f3; and the divergence is coordinate independent 
since (T MI/ )cM,ren(T/ji/)cM,ren als o diverge s there. Though the conformal Minkowski vacuum is a good vacuum for 
simply connected de Sitter space 1 158 15E], it (in the adapted version) is not self-consistent for the multiply connected 
de Sitter space dS/B. (That is, it does not solve the semiclassical Einstein's equations. 



In the case of an eternal Schwarzschild black hole, there are the Boulware vacuum [ 168 1 and the Hartle-Hawking 
vacuum [ 169 1 . The globally defined Hartle-Hawking vacuum bears essentially the same relationship to the Boulware 
vacuum as the Minkowski vacuum does to the Rindler vacuum [17C ] . For the Boulware vacuum, the energy- momentum 
tensor diverges at the event horizon of the Schwarzschild black hole, which means that this state is not a good 
vacuum for the Schwarzschild black hole because, when one inserts this energy-momentum tensor back into Einstein's 
equations, the back-reaction will seriously alter the Schwarzschild geometry near the event horizon. For the Hartle- 
Hawking vacuum, howev er, t he energy-momentum tensor is finite everywhere and a static observer outside the horizon 
sees Hawking radiation [171]. People usually regard the Hartle-Hawking vacuum as the reasonable vacuum state for 
an eternal Schwarzschild black hole because, when its e nergy-momentum tensor is fed back into Einstein's equations, 
the Schwarzschild geometry is only altered slightly [172]. Therefore, in the case of Misner space, Li and Gott [3(j have 
tried to find a vacuum which is also self-consistent and found that the adapted Rindler vacuum is such a vacuum if 
b = 2ir. 
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Here we also try to find a self-consistent vacuum for our multiply connected de Sitter space. Let us consider the 
adapted conformal Rindler vacuum in dS/B. The energy-momentum tensor of a conformally coupled scalar field 
in the adapted Rindler vacuum in Misner space is given by Eq. (32). Inserting Eq. J32] ) into Eq. ( |S8| ) and using 
Eqs. (^9|-^), we obtain the energy-momentum tensor for the adapted conformal Rindler vacuum of a conformally 
coupled scalar field in our multiply connected de Sitter space 



U40n 2 r 4 (1 - r 2 /r 2 
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9607r 2 r^ ' 



(95) 



where the coordinate system is the static coordinate system (t,r,9,(f>). Similarly, this result can also be analytically 
extended to the whole dS/B, though the static coordinates only cover region 1Z. We find that, if 



/3-2^r , 



(96) 



this energy-momentum tensor is regular on the whole space. [Eq. (96) corresponds to b — 2ir via Eq. (p3[)] . Otherwise 
both (T /J y )cR,ron and (T A " y )cR,rcn(7) t ^)cR,rcn diverge as the Cauchy horizon is approached. For the case /3 = 27rro, 
the energy-momentum tensor is 



(Tab) 



CR,ron 



9607rV(5 



gab. 



(97) 



which is the same as the energy-momentum tensor for the conformal Minkowski vacuum in the simply connected 
de Sitter space. 

The Euclidean section of our multiply connected de Sitter space is a four-sphere S 4 embedded in a five dimensional 
flat Euclidean space with those points related by an azimuthal rotation with angle (3/r being identified. There are 
conical singularities unless (3/tq — 2n. This may be regarded as a geometrical explanation of the self-consistent 
condition in (|96|). 

Similarly, our multiply connected de Sitter space solves the semiclassical Einstein's equations with a cosmological 
constant A and the energy-momentum tensor in Eq. ( |97| ) (and thus it is self-consistent) if 7q 
(if A = 0, we have the two solutions r 2 = 1/3607T and r = 00 p9[). 
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D. Particle Detectors in the Multiply Connected de Sitter Space 



It is well known that in the simply connected de Sitter space, an inertial particle detector perceives thermal radiation 
with the Gibbons-Hawking temperature [Eq. (|5(])] if the conformally coupled scalar field is in the conformal Minkowski 
vacuum |100 152]. Now we want to find what a particle detector perceives in the adapted conformal Rindler vacuum 
in our multiply connected de Sitter space. 

The response function of the particle detector is still given by Eq. ( |l7| ) . The Wightman function is obtained from 
the corresponding Hadamard function by Eq. (|5^). The Hadamard function for the conformally coupled scalar field 
in multiply connected de Sitter space is related to that in Misner space via Eq. ( p4] ) [with G(X, X') replaced by 
G^(X,X')]. The Hadamard function for the adapted Rindler vacuum in Misner space is given by Eq. (pl|). Inserting 

Eq. (|l|) [as into Eq. (§J) and using Eqs. ( |79|j83| ), we obtain the Hadamard function for the adapted conformal 
Rindler vacuum of the conformally coupled scalar field in our multiply connected de Sitter space 



G, 



sinl^y^l - r 2 /rl){\ - r' 2 /r%) [-(t - t' + n/3) 2 + r 2 ~f 2 



where X — (t, r, 9, <fi), X' = (t' , r', 9' , cj)'), and 7 is written in (t, r, 

1 



cosh 7 = 



^/(l-r 2 /r 2 )(l~r> 2 /r 2 ) 



1 1 - ^- [cos 9 cos 9' + sin 9 sin 9' cos(0 - <j)')} j 



(98) 



(99) 



The Wightman function is obtained from Eq. ( |98| ) via Eq. (|5C|). The Hadamard function given by Eq. (|98j) and the 
Wightman function obtained from that are defined in region 1Z in the multiply connected de Sitter space, but they 
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can be analytically extended to region T via the transformation in Eq. (169). However, it should be noted that as we 



make the continuation from TZ to T, \j (1 — r 2 /tq)(1 — r' /r%) should be continued to be — J (i 2 jr\ — l)(i' 2 jr 1 — 1) 



instead of +y (t 2 jr\ — l)(i' 2 /rQ — 1). This is because if we take \J\ — z 2 = i\j 1 z 2 — 1, we should also take yl — z' 2 = 
i\J z' 2 — 1 (instead of —is/ z' 2 — 1) [z and z' should be continued along the same path), thus \ (1 — z 2 )(z' 2 — 1) = 



1 z 2 — z' 2 — 1) = -J (z 2 — l)(z' 2 — 1). Using similar transformations, the results can also be continued to 

regions V and C (we do not write them out because we do not use them here). 

We consider particle detectors moving along three kinds of worldlines in our multiply connected de Sitter space: 
1. A particle detector moving along a worldline with r,9,(j) — constants in TZ. In such a case, on the worldline of 

the particle detector, 7 is zero and the Hadamard function is reduced to 

1 °° 1 

CR[ ' ' 2^ 2 (l-r 2 /rl) n ^ oo (t-t' + n(i) 2 
1 00 1 

- -5? £ 7 — ( 10 °) 



n=-oo (t — t' -\- nfi\J\ — r 2 /r^j 



where r — tJ 1 — r 2 /r 2 is the proper time of the particle detector. The corresponding Wightman function obtained 
from Eq. (50) is 



'crA 

n— — 00 



oiaM-si E -— — =r~Y- (101) 

It — r + npy 1 — r z /Tq — ze I 

where e is an infinitesimal positive real number. Inserting it into Eq. (|l7|), obviously the integration over At = t — t' 
is zero since all poles of the integrand are in the upper-half plane of complex At while the integration contour is 
closed in the lower-half plane. Therefore the response function •F(A-E') is zero and no particles are detected. All of 
these worldlines are accelerated, except for the one at r = 0. 

2. A particle detector moving along a geodesic with l,6,4> — constant in region T . In this region the Hadamard 
function is 



G U(x,x>) = --, J2 r ~ ~ 7 ' (102) 

-" n=-oosinh 7v /(i 2 /r 2 -l)(t' 2 /r 2 -l) + n(3) 2 + r 2 7 2 ] 



where 7 is given by 

cosh 7 = — = { — 1 + ~r [cos^cos^' + sin sin 6' cos(<^> — </>')]} 

V(?/fg-l)(? a /rg-l) 1 r ° J 



(103) 



[Eq. ( J102| ) and Eq. (103) are obtained from Eq. ( |98| ) and Eq. ( |99| ) via the transformation in Eq. ( |69| ) respectively.] On 
the worldline of the particle detector, the Hadamard function is reduced to 

00 

G ( cklt') = -— £ , 7 = . (104) 

2n n=-oosinh 7v /(i 2 /r 2 -l)(^ 2 /r 2 -l) {-n 2 /3 2 + r 2 7 2 ) 

and cosh 7 is reduced to 

cosh7 = -= ^Q" 1 = . (105) 



^{t 2 /r 2 -l){t' 2 /r 2 -l) 



Using the proper time r defined by Eq. (71), on the worldline of the particle detector cosh 7 and G^p can be written 
as 

cosh 2T + cosh At - 2 , 

cosh7 = — — — , (106) 

' cosh IT- cosh At ' v ; 
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and 



£ 3 



7 



sinh7(cosh2T — cosh Ar)(n 2 6 2 — 7 2 ) 



(107) 



where r > 0, r' > 0, Ar = (r — r')/ro, T = (r + r')/2ro, and 6 = /3/ro- The Wightman function is equal to one half 
of the Hadamard function with Ar replaced by Ar — ie [Eq. (pG)]. Thus the response function is 



(108) 



where 
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r o Jo 



dT 



dATe~ lAEroAT x 
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7 



sinh 7 (cosh2T— cosh Ar)(n 2 6 2 — 7 2 ) 



(109) 
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Now we consider the poles in the complex Ar plane of the integrand in the integral of T n (AE). The poles are given 
by the equation 



7 



±nb. 



(110) 



(It is easy to check that cosh2T = coshAr docs not give any poles.) From Eq. (|110| ) and Eq. ( |106| ), we have (we 
neglect the term ie, and at the end of the calculation we return it back to the expressions) 



Solutions to Eq. (Ill) are 



cosh 2T + cosh Ar - 2 = cosh nb (cosh 2T - cosh Ar) 



Ar = Ar n + i2mir = Ar„ 



where 



Ar„ = ±Arccosh 



(coshnfr- l)cosh2T + 2 
coshn6 + 1 



±2Arcsinh sinhTtanh 
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(111) 



(112) 



(113) 



where Arccoshz is the prin cipa l value of arccoshz, and here it is real (similarly for Arcsinhz). We ne ed to ch eck if 
all Ar„ m are roots of Eq. ( |llO|) , because the number of roots might increase as we go from Eq. (110) to Eq. ( 111 ). 
[E.g., for any integer m, x m = ±2 + imn solv es the equation cosh(2x) = cosh 4; but, only xq = +2 solves the equ ation 
2x — 4.] Ar„ is obviously a root of Eq. ( |110| ). The question is: as Ar goes from Ar„ to Ar„ + z2to7t, does Eq. ( |106| ) 
give the same 7 which is a real value [±nb; see Eq. (110)]? (Remember that arccosh z is a multi- valued complex 
function.) To answer this question, let Ar = Ar„ + i9 (where 9 is real). Then from Eq. (106) we have 



7 = arccosh 



cosh 2T + cosh Ar - 2 
cosh 2T - cosh Ar 



sinhT + sinh^ 
sinh T — sinh 



sinh T + sinh ALh. cos | + j cosh ^La. s i n £ 
= In - 



sinh T — sinh -^p- cos | 



i cosh s i n £ 



Z2 



where we have used arccoshz = ln(z + y/z 2 — 1). The real components of z\ and z-i are respectively 

(z\ ) = sinh T + sinh T " cos ^ , 



(114) 



(115) 



sinh T — sinh • 



Ar„ 



— — cos - . 

2 2 



By Eq. (113), we find that ^t(zi) and 5ft(z2) are always positive for any real 9. 
Ar n to Ar„ + i2mir, the arguments (the argument of a complex number z = \z\e' 



(116) 

This means that as Ar goes from 
^ is a) of z\ and z 2 do not change, 
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neither does the argument of z\/zi- The value of 7 remains in the same branch of In z as 9 var ies. Thus, for all 
Ar„ m = At„ + i2rmr, we have 7 = ±nb and Eq. (|llO| ) is satisfied. Therefore all Ar„ m in Eq. (112) are poles. 
The residues of the integrand in (109) at poles Ar„ m are (here ie is returned to the expressions) 

Res(Ar = i2rmr + ie, n = 0) = ^J e 2 ™-A£r 0j (n7) 

^ g2rrnrAErQ— iAEr Ar n 

Res(Ar = Ar„ + i2mn + ie, n ^ 0) = : — ; . . . — . (118) 

v ' ^ ; 4?r 2 (coshn6+ l)sinhAr„ v ; 

Then by the residue theorem (the contour for the integral is closed in the lower-half plane of complex At) we have 

dT r AE 



dT 2tt e 2 * AEr ° - 1 ' 

and 

dT n ^ sin(A£;ro|Ar n |) 



dT 7r(coshn6 + 1) sinh |At„| e 27TAEr ° - 1 ' 



(119) 



(120) 



The sin( AEr \ Ar n \ ) factor in Eq. (120) indicates that the n^O terms' contribution can be both positive (absorption 
by the detector) and negative (emission from the detector). We see that the contribution of the n = term is just the 
Hawking radiation with the Gibbons-Hawking temperature Tq_h = 1/2tiTo in the simply connected de Sitter space. 
The contribution of the n^O terms is a kind of "grey-body" Hawking radiation: the temperature is Tq-h, but its 
density or flux decreases as the universe expands (|Ar„| increases as the universe expands). The sum of all n =/= 
contributions is 

^ dfn = 1 1 \ - sin(A£r |AT ra |) 

^ dT 7T 2 e 2 "» AB - 1 ^ (coshn&+ l)smh|ArJ' 

In the case of b = 2ir (the self-consistent case), we have coshnfr ~ exp(|n|6)/2 > 1 (n ^ 0) and thus Ar„ ~ ±2T. 
Then 

V dT n _ 1 A sm(2AEr T) 

^ dT ~2ne 2 ^ AE -l sinh2T ' [ ' 

n=£0 

where A = 4^^° =1 (cosh2n7r + 1) _1 ~ 0.015. As T — > 00, the contribution of all n ^ terms decreases exponentially 
to zero. Thus, at events far from the Cauchy horizon in T , the particle detector perceives pure Hawking radiation 
given by the n = term. As T — > (near the Cauchy horizon), we have 

^dT n _Ar AE 

^ dT ~ 2tt e 2 * r ° AE - 1 ' [ ' 

This is a "grey-body" Hawking radiation with A ~ 1.5%. Near the Cauchy horizon the total radiation is the sum 
of a pure Hawking radiation (given by the n = term) and a "grey-body" Hawking radiation (given by all n =/= 
terms). The total intensity of the radiation near the Cauchy horizon is a factor of ~ 101.5% that of regular Hawking 
radiation, but its spectrum is the same as the usual Hawking radiation. 

3. A particle detector moving along a co-moving worldline in the steady-state coordinate system. Suppose the 
detector moves along the geodesic p,0,<f> = constants (such a worldline is a timelike geodesic passing through 1Z and 
into T) where p = (x 2 + y 2 + z 2 ) 1 / 2 and the proper time r are related to the static radius r by 

r = -r p/f}= pe T/ro . (124) 

The Cauchy horizon is at r = ro, or p — —rj = r$e~ T / Va . On the worldline of the detector the Hadamard function is 

1 00 1 

g£1(T,At) = ^—^ V k, (125) 

2vr 2 r 2 2LsinhAl n £- oo72 _^ +n ^ 2 

where At = (t — t')/tq, T = (t + r')/2ro, L = pe T /ro = r(T)/ro, 7 is given by 
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cosh 7 



and t — t' is related to T and At by 



cosh ■ 



t-t 1 
ro 



1-L 2 
\7l + L 4 -2L 2 cosh At' 



cosh At — L 2 
vT + L A - 2L 2 cosh At' 



(126) 



(127) 



By analytical continuation, Eqs. (125 127) hold in the whole region covered by the steady-state coordinates in de Sitter 
space. The Wightman function G + is equal to one half of with At replaced by At — it [Eq. (|50|)1. The response 
function is T{ AE) = E^L-oo-^n {AE) where 



1 pOC pOO 

T n (AE) = — dT dAre- lAEroAT 

,/ — oc J — oo 



2L sinh ^ 



7 Z 



(128) 



At-»At 



The poles of the integrand in the complex-AT plane are given by 

t-t' 



ro 



nb = ±7. 



This together with Eq. (126) and Eq. (127) leads to 



The roots of Eq. fll30|) are 



where 



(cosh At — L 2 ) cosh nb + sinh At sinh nb = 1 — L 2 . 



At = Ar±+i2?nn = At± to , 



At± = In 1 + 2/i 2 ± 2fiy/l + y? 



nb, 



(129) 



(130) 



(131) 



(132) 



where [i = sinh(n6/2). By carefully checking At^ in Eq. ( h31h , a s we did in case 2, we find that: (1) For L < 1 
(or pe T < r , i.e., in region 1Z), only Ar^g = At^ solve Eq. Q129D ; (2) for L > 1 (or pe T > r , i.e., in region T), 
only At ? +„ = At+ + i2mn solve Eq. ( |129| ) . (Here it is as sum ed that b > In 2 and t he se lf-consistent case with 6 = 2j r 
obviously satisfies this condition.) All other At's in Eq. ( |l3l| ) are not roots of Eq. (129), though they solve Eq. ( 130 ). 
Therefore the poles are (where ie is returned) 



At = 



Att 



ie, 



in 1Z: 



At+ + i2mn + ie, in T. 



Obviously in region 1Z all poles are in the upper-half plane of complex At. Therefore 

dT 



dT 



= 0, 



(133) 



(134) 



when the particle detector is in region 1Z. So the particle detector sees nothing while it is in region 7Z. 

In region T , only the poles with m < are in the lower-half plane of complex At. The residues of the integrand 
at poles A+ m + ie are 



Res(AT = i2mir + it, n = 0) = ^e 2 ' MA& », 



(135) 



Res(AT = Ar+ + i2m-K + it,n^ 0) 



167r 2 Lsinh^ 



ai(l + L 4 - 2L 2 coshAT+)e 2m * AEr °- iAEr ° AT i 
a\L(L 2 - l)cosh^i + (a 2 + n6) (L 2 cosh At ? | - 1) 



(136) 
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where a\ = Arccosh , L 1 and ai — Arccosh . L coshAr n gy ^he residue theorem, we have that 

V1 + L 4 -2L 2 cosh At+ Vl+L 4 -2L 2 coshA-r+ 



dJ-o/dT has the same value as that in Eq. ( |119| ), which represents Hawking radiation with the Gibbons-Hawking 
temperature; the contribution of all n ^ terms (note that A+ = — A+„) is 

d ^ 1 ^ sin(ASr Ar+) ^ 



a 1 (l + L 4 -2i 2 coshAr+) 



aiL(L 2 - 1) cosh - (a 2 + n6)(L 2 cosh At^ - 1) 



which represents a "grey-body" Hawking radiation. As T — ► oo (or L — * oo), ^ ^ tl ^ ^ exponentially drops to zero; 
therefore, at events far from the Cauchy horizon in the particle detector only perceives pure Hawking radiation 
(the same as that in case 2). As L — > 1 (approaching the Cauchy horizon), we also have gyXm^o^ 7 " ~ * 0- Thus 
as the Cauchy horizon is approached from the side of region the particle detector comoving in the steady-state 
coordinate system perceives pure Hawking radiation with Gibbons-Hawking temperature. 

From the above discussion, we find that in our multiply connected de Sitter space with the adapted Rindler vacuum, 
region 1Z is cold (where the temperature is zero) but region T is hot (where the temperature is Tq-h)- Similarly, 
region L is cold but V is hot, the above results can be easily extended to these regions. This gives rise to an arrow 
of increasing entropy, from a cold region to a hot region (Fig. |^). 



E. Classical Stability of the Cauchy Horizon and the Arrow of Time 

In classical electromagnetic theory, it is well known that both the retarded potential 4> rc t and the advanced potential 
0adv (and any part-retarded- and-part-advanced potential a<p le t + &</> a dv with a + b = 1) are solutions of Maxwell's 
equations. But from our experience, we know that all the electromagnetic perturbations we see are propagated only 
by the retarded potential. (For example, if at some time and some place, a light signal is emitted, it can only be 
received by a receiver at another place sometime later) . This indicates that there is an arrow of time in the solutions 
of Maxwell's equations, though Maxwell's equations themselves are time-symmetric. This arrow of time is sometimes 
called the electromagnetic arrow of time, or the causal arrow of time. How this arrow of time arises is a mystery. 



Many people have tried to solve this pr oblem b y attributing it to a boundary condition of the Universe ]60, 173-175 



(for revi ew of the arrows of time, see [ 176, 177| ) . In this subsection we argue that the principle of self-consistency 
[178,179| naturally gives rise to an arrow of time in our multiply connected de Sitter space. 

First let us consider the arrow of time in Misner space. Suppose at an event E in region F in Misner space [by 
boost and translation, assume we have moved E to (t = to,x = 0, y = 0, z — 0)], a spherical pulse of electromagnetic 
wave is created. If the potential is retarded [here "retarded" and "advanced" are defined relative to the direction of 
(d/dt) a (t is the time coordinate in the global Cartesian coordinates of the covering space — Minkowski space)], the 
pulse will propagate in the future direction as a light cone originating from E. At any point on the light cone, the 
energy-momentum tensor of the wave is 

T ab = fik a k b , (138) 

where fj, = fi(t) is a scalar function and k a = k°(d/dt) a + fc 1 (9/9x) a + k 2 (d/dy) a + k 3 (d/dz) a is a null vector tangent 
to the light cone, and the energy density measured by an observer with four-velocity vector (d/dt) a (whose ordinary 
three-velocity is zero) is 

'- r -(s)'(s)'-^ )1 - (139) 

(Thus fi measures the energy density of the electromagnetic wave.) By Einstein's equations, the back-reaction of 
Tab on R and R a bR ab (where R a b is the Ricci tensor and R = R a a is the Ricci scalar curvature) is SR ~ T Q °, 
6(Ra b R ab ) ~ T ab T ab . The Riemann tensor can be decomposed as R a bcd = C a bcd + Qabcd, where C a bcd is the Weyl 
tensor and Q a bcd is constructed entirely from the Ricci tensor 

Qabcd = g a [ C Rd]b ~~ 9b[ C Rd]a - -^R9a[c9d]b, (140) 



where square brackets denote antisymmetrization [ ISO | . The Weyl tensor describes the part of the curvature that is 



due to pure gravitational field, whereas the Ricci tensor describes the part that, according to Einstein's equations, is 
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directly due to the energy- momentum tensor of matter ]59| , Therefore, in some sense, the values of T a a and T ab T ab 
determine the influence of matter fields on the stability of the background spacetime. An infinite T a a or T ab T ab 
implies that the spacetime is unstable against this perturbation and a singularity may form; on the other hand, if T a a 
and T ab T ab are finite, the spacetime may be stable against this perturbation. Self-consistent solutions should require 
that T a a and T ab T ab do not blow up. If they did, the starting geometry — on the basis of which T a a and T ab T ab 
were calculated — would be greatly perturbed and the T a a and T ab T ab calculation itself would be invalid, and thus it 
would not be a self-cons isten t solution. For electromagnetic fields we always have T a a = 0, so we need only consider 
T ab T ab . For T ab in Eq. (Egf), we also have 



a b 



o. 



(141) 



Thus significant perturbations (indicated by a non-vanishing T ab T ab ) can only occur when the light cone "collides" 
with its images under the boost transformation. At any point p on the intersection of the light cone L and its n-th 
image L n (suppose n > 0), the energy- momentum tensor is 



T ab = uk a k b + jlk a k b 



(142) 



where fc° is the null vector tangent to the light cone L at p, k a is the null vector tangent to the light co ne L n at p; p 
measures the energy density in light cone L, p measures the energy density in light cone L n . From Eq. (142) we have 



T ab T ab = [2pp(k a k a ) 2 ] p , 



(143) 



the index p denotes that the quantity is evaluated at the point p. 

Since the point p on L n is obtained from some point p 1 on L by boost transformation, p and p 1 must have the same 
timelike separation from the origin (t = 0, x = 0, y — 0, z — 0) (remember that p is on the intersection of L and L n , 
see Fig. ^|a). If we take the k a at p being transported from the k' a at p', we have A* p gl„ = Vp'eL- Because the light 
cone L is spherically symmetric, we have t p =t p >. Therefore we have /J, p >eL = l^peL and at p we have p, = /i. Under 
the boost transformation B, we have 



(k a ) P = B[(k' a ) p/ ] = fc' C 



k' 1 



coshn& 



d_ 

d.r 



coshnfr f — 
f sinh nb 



sinh nb 
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dx 



k' 2 



k" 
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where {k' a ) p , = k'°(d/dt) a + k'\d/dx) a + k' 2 (d/dy) a + k' 3 (d/dz) a . Due to the spherical symmetry, we have k'° = k° 
Define {r,6,(j)) by x = rcosO, y = rsin#cos0, and z = rsin^sint/). Then we have r' — r, 9' = it — 0, <f> 
means "at p'"), and 



k 1 = k a cos 8, k 2 = k° sin 9 cos < 



k° sin 9 sin < 



("/" 

(145) 



and 



Then 



and 



k' 1 = k'°co$e' = -fc cos6» = -k 1 , k' 2 = k'°smd'co8(f>' = fc°sinflcos^= k 2 
k' 3 =k'° sin 9' sin </>' = k° sin 9 sin = fc 3 . 



(k a k a ) p = (fc ) 2 [- (1 + cos 2 9) cosh nb + 2 cos 9 sinh nb + sin 2 9} , 



T ab T ab = 2p 2 {t p )[-{l 



cos 2 9) cosh nb + 2 cos 9 sinh nb + sin 2 01 2 . 



It is easy to find that T T ab reaches a maximum at 9 — and 



(T ab T ob ) r 



8p 2 (t p )e 



-2nb 



(146) 
(147) 
(148) 
(149) 



where p(t p ) is the energy density from L as measured in a frame at event p with ordinary velocity v x — v y = v z = 0. 
(T oh T ofc ) max is always finite [less than 8p 2 (t p )] since n is positive. If n < we have (T afc T afc ) max = 8p 2 (t p )e 2nb < 
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8p 2 (t p ). So if we have a retarded potential in region F, even considering the infinite number of images, T ab T ab is 
always finite. 

If the potential is advanced however, the pulse wave will propagate backward in the past direction as a light cone 
originating from E. And, within a finite time, it will hit the Cauchy horizon. By an analysis similar to the above 
arguments, we find that in this case 

T ab T ab = 2p 2 (t p )[(l + cos 2 9) coshnb + 2 cos 6 sinh nb - sin 2 9} 2 , (150) 

which reaches a maximum at 9 = and 

(T ab T ab ) max = 8p 2 (t p )e 2 ^ b . (151) 

Since p(t p ) is finite (the past light cone from E at 8 — hits the Cauchy horizon in a finite affine distance), thus 
(T T a b) max — > oo as n — » ±oo. As n — ► ±oo, L and L n collide at the Cauchy horizon [as n — * ±oo the point p(9 = 0) 
approaches the Cauchy horizon] (see Fig. ^p). Thus {T ab T ab ) max diverges as the Cauchy horizon is approached and 
the Cauchy horizon may be destroyed. Therefore the advanced potential is not self-consistent in region F of Misner 
space. It is easy to see that any part-retarded-and-part-advanced potential is also not self-consistent in F. The only 
self-consistent potential in region F is the retarded potential. 

Similarly, in region P the only self-consistent potential is the advanced potential (see Fig. ^]a). [Note that here 
"advanced" and "retarded" are defined relative to the global time direction in Minkowski spacetime (the covering 
space) . An observer in P will regard it as "retarded" relative to his own time direction.] 

In region R, by boost and translation, we can always move the event E (where a spherical pulse of electromagnetic 
waves is emitted) to (i = 0, x = x$, y = 0, z = 0). Either pure retarded or pure advanced potentials are self-consistent 
in this region because the light cone never "collides" with the images of itself and thus we always have T ab T ab = 
(see Fig. gb). But, for a part-retarded-and-part-advanced potential, the retarded light cone (L + ) propagates forward 
while the advanced light cone (L~) propagates backward, both originating from E. The forward part of the light cone 
will collide with images of the backward part of the light cone and vice versa (see Fig. ^d) . We find that at a point p 
on the intersection of L + and L~ (or L~ and L+) 

T ab T ab = 2p(t)p{-t)[(l + cos 2 8) cosh nb- 2 cos 6» sinh nfe + sin 2 6»] 2 , (152) 

where p(t) is the energy density from L + observed in a frame on L + with time coordinate t and with ordinary velocity 
v x — Vy — v z — and p{—t) is the energy density from L~ seen in a frame on L~ with time coordinate —t and with 
ordinary velocity v x — v y = v z — 0. T ab T ab reaches a maximum at 9 = ir, and 

(T ab T ab ) max = 8p(t)p(~t)e 2 ^ b , (153) 

where t is the global time coordinate in the covering Minkowski space. As p approaches the Cauchy horizon, where 
n — > ±oo, p(t) and p{—t) are both finite, since in the 9 = tt direction the future and past light cones of E both 
hit the Cauchy horizon in a finite affine distance. Thus {T ab T ab ) max — ► oo as p approaches the Cauchy horizon 
(where n — > ±oo). Therefore in region R both the retarded and the advanced potential are self-consistent, but the 
part-retarded-and-part-advanced potential is not self-consistent. This conclusion also holds for region L. Furthermore, 
there must be a correlation between time arrows in region L and region R: if we choose the retarded potential in R, we 
must choose the advanced potential in L (see Fig. ^b); if we choose the advanced potential in R, we must choose the 
retarded potential in L. Otherwise the collision of light cones from R and light cones from L will destroy the Cauchy 
horizon. 

As another treatment for perturbations in Misner space, consider that at an event E in region F two photons are 



created [181] [we choose E to be at (t = to, x = 0, y = 0, z = 0) as before]. One photon runs to the right along the 
+x direction, the other photon runs to the left along the ~x direction. They have the same frequency (thus the same 
energy). The tangent vector of the null geodesic of the right-moving photon is chosen to be r k a — f (^;) a = (gf~) a j 
where A r is an affine parameter of the geodesic, q is a constant and u — t + x, v — t — x. The tangent vector of 
the null geodesic of the left-moving photon is chosen to be ik a = ^{-§^) a = where A; is an affine parameter 

of that geodesic. The null vectors r k a and ik a are invariant under boost transformations. At any point where a 
photon with null wave- vector k a is passing by, the frequency of the photon measured in a frame of reference passing 
by the same point with the four-velocity v a is oj = —k a v a . If v a = (d/dt) a (i.e., the frame of reference has ordinary 
three- velocity v x = v y = v z = 0) and k a = r k a or ;fc a , we have Lu r — u>i — q/2to = luq (thus q measures the frequency 
of the photon). At any point where the n-th image of the right-moving (left- moving) photon is passing by, using 
the boost transformation we can always find a frame of reference in which the frequency of the photon is u>q. But 
at a point p where the right-moving (left-moving) photon passes the n-th image of the left-moving (right-moving) 
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photon, we cannot find a frame of reference such that the two "colliding" photons both have frequency ujq. In such a 
case we should analyze it in the center-of-momentum frame. The four-velocity of the center-of-momentum frame is 
v a = 7( r fc a + ik a ) where j 2 — — [{ r k a + ik a )( r k a + j/ca)]" 1 = uv/q 2 — ff/q 2 where rj = (t 2 — x 2 ) 1 ^ 2 is the proper 
time separation of p from the origin (t = 0, x = 0, y = 0, z = 0). Therefore the total energy of the two oppositely 
directed photons in the center-of-momentum frame is 

1 It 

£ = wi + uj 2 = - r k a v a - ik a v a = - = —lo . (154) 

7 V 

(For all other frames the total energy would be greater.) If the potential is retarded, so photons move in the future 
direction, all points where photons and their images "collide" are in the future of the hypersurface t 2 — x 2 = t 2 y 
Therefore we have rj > fjo = to and £ < 2wo, so the total energy in the center-of-momentum frame is always bounded. 
But, if the potential is advanced, photons move in the past direction; thus all points where photons and oppositely 
directed image photons "collide" are in the past of the hypersurface t 2 — x 2 — tg. In particular, the right-moving 
(left-moving) photon collides with the oo-th (— oo-th) image of the left-moving (right-moving) photon at the Cauchy 
horizon, where rj = and thus £ — > oo. Thus, the Cauchy horizon may be destroyed by these photon pairs. Therefore 
in agreement with our earlier argument, the advanced potential is not self-consistent in region F. The retarded potential 
is self-consistent in region F. Similarly, the advanced potential is self-consistent in region P. In region R and region 
L, both the retarded potential and the advanced potential are self-consistent, because the photons and their images 
will not collide with each other and at any point a photon is passing by we can always find a frame for whom the 
frequency of this photon is luq. And, the potentials in region R and region L must be correlated in the following 
way: If the potential in R is retarded, the potential in L must be advanced; if the potential in R is advanced, the 
potential in L must be retarded (we would call them "anti-correlated"). Otherwise the photons from L and photons 
from R passing in opposite directions would be measured to have infinite energy in center-of-momentum frames as the 
Cauchy horizon is approached and this may similarly destroy the Cauchy horizon. These conclusions are consistent 
with those obtained from the analysis of the perturbation of a pulse wave discussed above. 

Our multiply connected de Sitter space is conformally relat ed to Misner space via Eqs. ( [78[j83] ). Because light 
cones and chronological relations are conformally invariant |l82[ (thus regions J 7 , V, 7Z, and C in multiply connected 
de Sitter s pace c orrespond respectively to regions F, P, R, and L in Mis ner space under the conformal map, as discussed 



in section KB ), Maxwell's equations are also conformally invariant [18C ,182 j, so it is easy to generalize the results 
from Misner space to our multiply connected de Sitter space. Under the conformal transfo rmat ion g ab — ► fl 2 g abl 
the energy-momentum tensor of the electromagnetic field is transformed as T a b — ► il~ 4 T a b |l80| . Thus T ab T ab is 
transformed as T ab T ab -> n- s T ab T ab . From the above discussion of T ab T ab in Misner space, we know that T ab T ab is 
zero everywhere except at the intersection of two light cones. Thus, in multiply connected de Sitters pace, T ab T ab is 
also zero everywhere except at the intersection of two light cones. At the intersection of two light cones in multiply 
connected de Sitter space, it is easy to show that the maximum value of T ab T ab is at the points with 8 — or 8 = ir on 
the intersection. From Eq. (^0|) and Eq. (|8^) we find that for 8 = or 8 = n, fl 2 is non-zero except at the points with 
8 = on the Cauchy horizon (where r = ro or t = ro). Also because fi 2 is finite everywhere on the Cauchy horizon (i.e. 
it is never infinite), we have that: (1) if T ab T ab diverges on the Cauchy horizon in Misner space, the corresponding 
T ab T ab also diverges on the Cauchy horizon in our multiply connected de Sitter space; (2) if T ab T ab is finite in some 
region (except at the Cauchy horizon) in Misner space, the corresponding T ab T ab is also finite in the corresponding 
region (not at the Cauchy horizon) in the multiply connected de Sitter space; (3) if T ab T ab is zero in some region 
(not a single point) in Misner space, the corresponding T ab T ab is also zero in the corresponding region in the multiply 
connected de Sitter space. Under the conformal transformation g ab — > £l 2 g ab , the affine parameter of a null geodesic 
is transformed as A — * A : dX/dX = Cfl 2 where C is a constant ]180| and thus the null vector k a — {d/dX) a is 
transformed as k a — > C^ 1 Qr 2 k a . Then 7 = [— ( r k a + ik a ){ r k a + ik a )]~ L l 2 is transformed as 7 — ■> Cfl^f and the total 
energy of the photon pairs in the center-of-momentum frame is transformed as £ — > C~ 1 fl~ 1 £ and the constant C _1 
can be absorbed into luq. Therefore, we can transplant the above results for Misner space directly to our multiply 
connected de Sitter space: In region T the only self-consistent potential is the retarded potential; in region V the only 
self- consistent potential is the advanced potential; in regions 1Z and C both the retarded potential and the advanced 
potential are self-consistent, but they must be anti-correlated (Fig. ^|). 

The Cauchy horizon separating a region with CTCs from that without closed causal curves is also called a 
chronology horizon || . A chronology horizon is called a future chronology horizon if the region with CTCs lies to the 
future of the region without closed causal curves; a chronology horizon is called a past chronology horizon if the region 
with CTCs is in the past of the region without closed causal curves. It is generally believed that a future chronology 
horizon is classically unstable unless there is some diverging effect near the horizon J9|,[50| . The argument says that a 
wave packet propagating in the future direction in this spacetime will pile up on the future chronology horizon and 
destroy the horizon due to the effect of the infinite blue-shift of the frequency (and thus the energy) seen by a timelike 



33 



observer near a closed null geodesic on the horizon [181,||. But if there is some diverging mechanism (like the diverging 



effect of a wormhole in a spacetime with CTCs constructed from a wormhole 0) near the horizon, the amplitude of 
the wave packet will decrease with time due to this mechanism, and this may cancel the effect of the blue-shift of the 
frequency, making the energy finite and thus rendering the future chronology horizon classically stable. Unfortunately, 
in our multiply connected de Sitter spacetime (as also in Misner space) there is no such diverging mechanism. A ligh t 



ray propagating in de Sitter space will focus rather than diverge. This can be seen from the focusing equation [ 183 1 



(155) 



where A is the cross-sectional area of the bundle of rays, A is the affine parameter along the central ray, the null 
vector k a is k a = (d/d\) a , and a is the magnitude of the shear of the rays. For de Sitter space we have R a bk a k b = 

h-gabk a k b ~ and thus we have d ^ < 0, so the ray will never diverge. (In fact this always holds if the spa cetim e 



satisfies cither the weak energy condition or the strong energy condition and it is called the focusing theorem [ 183 1 .) 
Hawking [pot has given a general proof along the above lines that any future chronology horizon is classically unstable 
unless light rays are diverging when they propagate near the chronology horizon. You could cause this instability by 
shaking an electron in the vicinity of the future chronology horizon. The retarded wave would then propagate to the 
future causing the instability. 

However, in Hawking's proof pp| , if we replace a future chronology horizon with a past chronology horizon, then 
the proof breaks down because, in such a case, a wave packet propagating toward the future near the past chronology 
horizon will suffer a red-shift instead of a blue-shift. Therefore a past chronology horizon, according to Hawking's 
argument, is classically stable in a world with retarded potentials. If the universe started with a region of CTCs, but 
there are no CTCs now, that early region of CTCs would be bounded to the future by a past chronology horizon, 
and that horizon would be classically stable in a world with retarded potentials — which is what we want. In our 
multiply connected de Sitter space, this is realized, since the arrow of time in region T is in the future direction and 
the arrow of time in region V is in the past direction [here "future" and "past" are defined globally by the direction 
of (d/dr) 0, , where r is the time coordinate in the global coordinate system (r, x, 9, (j>) of the de Sitter covering space]. 
T and TZ can have retarded potentials, while V and C have advanced potentials, as we have noted. In this case the 
Cauchy horizons separating T from TZ and V from C are classically stable, as indicated by our detailed study of 
T ab T a t, as these Cauchy horizons are approached. What about the Cauchy horizons separating V from TZ and T from 
CI In region V ', the potentials are advanced, so Hawking's instability does not arise as one approaches the Cauchy 
horizon separating it from TZ. In region TZ, the potentials are retarded, so by Hawking's argument, one might think 
that there would be an instability as the Cauchy horizon separating 1Z from V is approached from the 1Z side. But, 
as we have shown, with retarded potentials in 1Z, T ab T a b does not diverge as the Cauchy horizon separating 1Z from 
V is approached from the 1Z side, indicating no instability. Why? Because one can always find frames where the 
passing photon energies are bounded as the Cauchy horizon is approached. Hawking's argument works only if one 
can pick a particular frame like the frame of a timelike observer crossing the Cauchy horizons and observe the blow 
up of the energy in that frame. (Thus Hawking's approach is observer-dependent, while our approach with T ab T a b is 
observer-independent.) Hawking's timelike observer would be killed by these photons. But, as we have shown, 7Z is in 
a pure vacuum state in our model, so there are no timelike observers in this region, and no preferred frame. If there 
were timelike particles of positive mass crossing from V to 1Z through the Cauchy horizon, we have shown (Li and 
Gott |p0|) that these would cause a classical instability; but there are none. There are, as we shall show in the next 
subsection, no real particles in regions C and TZ (because these are vacuum states) and no real particles in region T 
and V until the vacuum state there decays by forming bubbles at a timelike separation |t| > To from the origin (tq 
will be given in the next subsection). Thus, there are no particles crossing the Cauchy horizons separating V from TZ 
and T form C Thus, there is no instability caused by particles crossing the Cauchy horizons; and since there are no 
timelike observers in region TZ to be hit by photons as the Cauchy horizon separating TZ from V is approached, there 
is no instability, as indicated by the fact T ab T a b does not blow up as that Cauchy horizon is approached. As indicated 
in Fig. ^b, region T + TZ is one causally connected region which can be pictured as partially bounded to the future by 
the future light cone of an event E' and bounded to the past by the future light cone of an event E; but E and E' are 
identified by the action of the boost, so these two light cones are identified, creating a periodic boundary condition 
for region T + TZ. As our treatment using T a bT ab with images indicates, retarded photons created in T + TZ cause no 
instability. Particles with timelike worldlines crossing the Cauchy horizons separating T + TZ from V + C would cause 
instability by crossing an infinite number of times between the future light cones of E and E', thus making an infinite 
number of passages through the region T + TZ (also V + C) shown in Fig. |]b. However, as we have shown, there 
should be no such particles with timelike worldlines crossing the Cauchy horizons separating T + TZ from V + C, and 
no photons crossing these horizons cither, since the potentials in T + TZ are retarded, while the potentials in V + L 
are advanced. Thus, we expect T + TZ and V + C to both be stable, and causally disconnected from each other. (See 
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further discussion in the next subsection). 

Thus, the principle of self-consistency ||l 78 ,1791 produces classical stability of the Cauchy horizons and naturally 



gives rise to an arrow of time in our model of the Universe. 

F. Bubble Formation in the Multiply Connected de Sitter Space 

From the above discussion we find that in the multiply connected de Sitter space region T and region V are causally 
independent in physics: the self-consistent potential in T is the retarded potential, while the self-consistent potential 
in V is the advanced potential, thus an event in T can never influence an event in V, and vice versa. T and V are 
physically disconnected though they are mathematically connected. If we choose the potential in TZ to be retarded, 
then the potential in C must be advanced. (Note that here "advanced" and "retarded" are defined relative to the 
global time direction in de Sitter space — the covering space of our multiply connected de Sitter space.) Then region 
T + TZ (including the Cauchy horizon separating T from TZ) forms a causal unit, and region V + C (including the 
Cauchy horizon separating V from C) forms another causal unit. (See Fig. 4b, where the two null surfaces partially 
bounding the grey T + TZ region to the past and future are identified. Similarly for the null surfaces partially bounding 
the V + £ region.) An event in J-+TZ and an event in V + £ are always causally independent in physics: they can never 
physically influence each other though they may be mathematically connected by some causal curves (null curves or 
timclike curves). Though T + TZ and V + C are connected in mathematics, they are disconnected in physics. They 
are separated by a Cauchy horizon. When we consider physics in JF + TZ, we can completely forget region V + C (and 
vice versa). Though in such a case the Cauchy horizon separating J-' + TZ from V + £ is a null spacetime boundary, we 
do not need any boundary condition on it because the topological multi-connectivity in T + TZ has already given rise 
to a periodic boundary condition (which is a kind of self-consistent boundary condition). (In Fig. 4b this is shown by 
the fact that the null curves partially bounding T + TZ to the past and future are identified.) This periodic boundary 
condition (the self-consistent condition) is sufficient to fix the solutions of the universe. For example, in our multiply 
connected de Sitter space model, the stability of the Cauchy horizon requires that the regions with CTCs (TZ and C) 
must be confined in the past and in these regions all quantum fields must be in vacuum states (as we have already 



remarked, the appearance of any real particles there seems to destroy the Cauchy horizon 30 1). This gives rise to an 
arrow of time and an arrow of entropy in this model. 

J-+TZ is a Hausdorff manifold with a null boundary, and thus J-+TZ is geodesically incomplete to the past. But, the 
geodesic incompleteness of T -\-TZ may not be important in physics because in the inflationary scenario all real particles 
are created during the reheating process after inflation within bubbles created in region T and these particles emit 
only retarded photons which never run off the spacetime because here the geodesic incompleteness takes place only in 
the past direction. On the other hand, we can smoothly extend T + TZ to V + C so that the total multiply connected 
de Si tter space dS/B is geodesically complete but at the price that it is not a manifold at a two-sphere (section 



IX A). This model describes two physically disconnected but mathematically connected universes. [The analogy 



between the causal structures in region T + TZ and region V + C might motivate us to identify antipodal points in our 



multiply connected de Sitter space, as we did for the simply connected de Sitter space (section VIII). The spacetime 
so obtained is a Hausdorff manifold everywhere. It is geodesically complete but not time orientable. For computing 
the energy-momentum tensor of vacuum polarization, we must take into account the images of antipodal points in 
addition to the images produced by the boost transformation. Further research is needed to find a self-consistent 
vacuum for this spacetime.] 

Now we consider formation of bubbles in T + TZ in multiply connected de Sitter space. [The results (and the 
arguments for T + TZ in the previous paragraph) also apply to region V + C, except that while in T + TZ bubbles 
expand in the future direction, in V + C they expand in the past direction; here "future" and "past" are defined with 
respect to (<9/<9r) a where r is the time coordinate in the global coordinates of de Sitter space.] Region TZ (for its 
fundamental cell see Fig. ||) which is multiply connected has a finite four- volume V\ = §7r&r 4 (here b = (3/ 'ro, f3 is the 
de Sitter boost parameter). If the probability of forming a bubble per volume r 4 in de Sitter space is e, then the total 
probability of forming a bubble in V\ is Pi = |7r6e. 

Region J- (its fundamental cell is shown in Fig. |J) has an infinite four-volume and thus there should be an infinite 
number of bubble universes formed [^9|,[)8| . The metric in region T is given by Eq. @ with 0<r<oo, <l < P, 
< 8 < 7r, and < <fi < 2ir (see Fig. 4a); it is multiply connected (periodic in / with period (3). In order that the 
inflation proceeds and the bubbles (which expand to the future — as expected with the retarded potential in region 
T) do not percolate, it is required that e < e cr where 5.8 x 1CP 9 < e cr < 0.24 B3]. Gott and Statler |^| showed 
that in order that we on earth today should not have witnessed another bubble colliding with ours within our past 
light cone (with 95% confidence) e must be less than 7.60 x 10" 4 for f2 = 0.1 (for Q = 0.4 Gott [§| found e < 0.01). 
In our multiply connected de Sitter space, for inflation to proceed, there should be the additional requirement that 
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bubbles do not collide with images of themselves (producing percolation) . A necessary condition for a bubble formed 
in T not to collide with itself is that from time r when the bubble forms to future infinity (r — > oo) a light signal 
moving along the I direction [where r and I are defined in Eq. (|6^) and Eq. fl7l])] propagates a co-moving distance 
less than f3/2, which leads to the condition that r > To = ro In ■ I n f a °t this is also a sufficient condition, 

which can be shown by the conformal mapping between region T in the multiply connected de Sitter space and region 
F-0 in Misner space defined by Eqs. ( pl[ - |83| ). If the collision of two light cones in F occurs beyond the hyperbola 
t 2 — x 2 — y 2 — z 2 = 1 (t > 0) in Misner space (i.e., in the region O), the corresponding two light cones (and thus 
the bubbles formed inside these light cones) in T will never collide because t 2 — x 2 — y 2 — z 2 = 1 in F corresponds 
to t — > 00 in T. It is easy to show that the condition for a light cone not to collide with its images within F-0 
is that e b (t 2 — x 2 ) — y 2 — z 2 > 1, where (t,x,y,z) is the event where the light cone originates. By Eq. ( |Sl| ) this 
condition corresponds to e b [(^) 2 — 1] > 1 + (^) 2 — 2^ cos^. Since t > ro and —1 < cos9 < 1, a sufficient condition 
is e b [(X) 2 - 1] > 1 + (X) 2 + 2±, i.e. e b (^ - 1) > X + 1 w hich is equivalent to r > r = r In Therefore all 

bubbles formed after the epoch To in T in the multiply connected de Sitter space will never collide with themselves. 
The < t < ro part of the fundamental cell in T has a finite four- volume Vn = Ui(cosh 3 ^ — 1). The total probability 

of forming a bubble in Vn is P\\ — |7r6e(cosh 3 p — 1). For b — 2tt we have To ~ 0.086ro, Vn ~ O.OIIV11, and thus 
Pii~0.011Pi. 

For the case of b = 2ir, in order that there be less than a 5% chance that a bubble forms in Vi (and thus less than 
0.05% chance in Vn), e should be less than 2 x 10~ 3 . This should be no problem because we expect that this tunneling 
probability e should be exponentially small. Thus it would not be surprising to find region 1Z and region T for epochs 
< t < To = 0.086ro clear of bubble formation events (and clear of real particles), which is all we require. 

Also note that there may be two epochs of inflation, one at the Planck scale caused by (T a f,) ren = — g a h/9607r 2 ro 
[Eq. (p7|)] which later decays in region T at t 3> to into an inflationary metastable state at the GUT scale produced 
by a potential V ((/>), which, still later, forms bubble universes. 



X. BABY UNIVERSE MODELS 



Inflationary universes can lead to the formation of baby universes in several different scenarios. If one of these 
baby universes simply turns out to be the original universe that one started out with, we have a multiply connected 
solution in many ways similar to our multiply connected de Sitter space. There would be a multiply connected region 
of CTCs bounded by a past Cauchy horizon which would be stable because of the self-consistency requirement as in 
the previous section, and this would also engender pure retarded potentials. Thus, in a wide class of scenarios, the 
epoch of CTCs would be long over by now, as we would be one of the many later-formed bubble universes. Also, 
the model might either be geodesically complete to the past or not. This might not be a problem in physics since 
we would in any case have a periodic boundary condition; and because with its pure retarded potentials, no causal 
signals could be propagated to the past in any case. There are several different baby universe scenarios — any one of 
which could accommodate our type of model. 

First, there is the Farhi, Guth, and Guven [3^] method of creation of baby universes in the lab. At late times in 
an open universe, for example, an advanced civilization might implode a mass (interestingly, it does not have to be 
a large mass — a few kilograms will do) with enough energy to drive it up to the GUT energy scale, whereupon it 
might settle into a metastable vacuum, creating a small spherical bubble of false vacuum with a V = A/8ir metastable 
vacuum inside. This could be done either by just driving the region up over the potential barrier, or by going close to 
the barrier and tunneling through. The inside of this vacuum bubble would contain a positive cosmological constant 
with a positive energy density and a negative pressure. This bubble could be created with an initial kinetic energy 
of expansion with the bubble wall moving outward. But the negative pressure would pull it inward, and it would 
eventually reach a point of maximum expansion (a classical turning point), after which it would start to collapse and 
would form a black hole. But occasionally, (probability P = 10~ 10 for typical GUT scales (32)) when it reaches its 
point of maximum expansi on it tunnels to a state of equal energy but a different geometry, like a doorknob, crossing 
the Einstein- Rosen bridge 184]. The "knob" itself would be the the interior of the bubble, containing the positive 
cosmological constant, and sitting in the metastable vacuum state with V = A/8tt. The "knob" consists of more 
than a hemisphere of an initially static S 3 closed de Sitter universe, where the bubble wall is a surface of constant 
"latitude" on this sphere. At the wall, the circumferential radius is thus decreasing as one moves outward toward the 
external spacetime. Just outside the wall is the Einstein-Rosen neck which reaches a minimum circumferential radius 
at r = 2M, and then the circumference increases to join the open external solution. This "doorknob" solution then 
evolves classically. The knob inflates to form a de Sitter space of eventually infinite size. It is connected to the original 
spacetime by the narrow Einstein-Rosen bridge. But an observer sitting at r — 2M in the Einstein-Rosen bridge will 
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shortly hit a singularity in the future, just as in the Schwarzschild solution. So the connection only lasts for a short 
time. The interior of the "knob" is hidden from an observer in the external spacetime by an event horizon at r = 2M . 



Eventually the black hole evaporates via Hawking radiation 185 1, leaving a flat external spacetime (actually part of 
an open Big Bang universe) with simply a coordinate singularity at r — as seen from outside. (See Fig. ^.) 

From the point of view of an observer sitting at the center of V = A/8ir bubble, he would see himself, just after the 
tunneling event, as sitting in a de Sitter space that was initially static but which starts to inflate. Centered on this 
observer's antipodal point in de Sitter space, he would see a bubble of ordinary V = vacuum surrounding a black 
hole of mass M. The observer sees his side of the Einstein- Rosen bridge and an event horizon at r = 2M which hides 
the external spacetime at late times from him. From the point of view of the de Sitter observer, the black hole also 
evaporates by Hawking radiation, eventually leaving an empty V = bubble in an ever-expanding de Sitter space. 
This infinitely expanding de Sitter space, which begins expanding at the tunneling event, is a perfect starting point 
(just like Vilenkin's tunneling universe) for making an infinite number of bubble universes, as this de Sitter space 
has a finite beginning and then expands forever. Now suppose one of these open bubble universes simply turns out 
to be the original open universe where that advanced civilization made the baby de Sitter universe in the first place 
(Fig. 0). Now the model is multiply connected, with no earliest event. There is a Cauchy horizon (C7i, see Fig. 0) 
separating the region of CTCs from the later region that does not contain them. This Cauchy horizon is generated by 
ingoing closed null geodesies that represent signals that could be sent toward the black hole, which then tunnel across 
the Euclidean tunneling section jumping across the Einstein-Rosen bridge and then continuing as ingoing signals to 
enter the de Sitter space and reach the open single bubble in the de Sitter space (that turns out to be the original 
bubble in which the tunneling event occurs). A retarded photon traveling around one of those closed null geodesies 
will be red-shifted more and more on each cycle, thus not causing an instability. Another novel effect is that although 
these null generators are converging just before the tunneling event, they are diverging just after the tunneling event, 
having jumped to the other side of the Einstein- Rosen bridge. Thus, converging rays are turned into diverging rays (as 
in the wormhole solution) during the tunneling event without violating the weak energy condition. These closed null 
geodesies need not be infinitely extendible in affine distance toward the past. It would seem that it can be arranged 
that the renormalized energy-momentum tensor does not blow up on this Cauchy horizon so that a self-consistent 
solution is possible. Using the method of images, note that the iV-th image is from N cycles around the multiply 
connected spacetime. The path connecting an observer to the N-th image will have to travel N times through the 
hot Big Bang phase which occurs in the open bubble after the false vacuum with V(<p < 4>o) = A/8ir dumps its false 
vacuum energy into thermal radiation as it falls off the plateau and reaches the true vacuum V{4> = <po) — 0. Thus, 
to reach the N-th image one has to pass through the hot optically thick thermal radiation of the hot Big Bang N 
times. And this will cause the contribution of the iV-th image to the renormalized energy-momentum tensor to be 
exponentially damped by a factor e~ Nr where r = nza t 1 (where n is the number density of target particles, z is 
the thickness of hot material, at is the total cross-section). Li p6| h as calculated the renormalized energy-momentum 
tensor of vacuum polarization with the effect of absorption. Li p6| has estimated the fluctuation of the metric of the 
background spacetime caused by vacuum polarization with absorption, which is a small number in most cases. If the 
absorption is caused by electron-positron pair production by a photon in a photon-electron collision, the maximum 
value of the metric fluctuation is (5g^ u ) msx ~ lp/(r c L), where lp is the Planck length, r c is the classical radius of 
electron, L is the spatial distance between the identified points in the frame of rest relative to the absorber |2(J . If we 
take L to be the Hubble radius at the recombination epoch (~ 10 23 cm), we have (5<7 Atl /) m ax ~ 10 -76 . Thus, we expect 
that the renormalized energy- momentum tensor will not blow up at the Cauchy horizon |2(f | , so that a self-consistent 
solution is possible. 

The tunneling event is shown as the epoch indicated by the dashed line in Fig. ^. During the tunneling event, the 
trajectory may be approximated as a classical space with four spacelike dimensions solving Einstein's equations, with 
the potential inverted, so that this Euclidean section bridges the gap between the two classical turning points. (In 
such a case, the concepts of CTCs and closed null curves should be generalized to contain a spacelike interval. Thus, 
there are neither closed null geodesies nor closed timelike geodesies with the traditional definitions. According to Li 
p3j, this kind of spacetime can be stable against vacuum polarization.) 

As Farhi, Guth, and Guven |32| note, the probability for forming such a universe is exponentially small, so an 
exponentially large number of trials would be required before an intelligent civilization would achieve this feat. If the 
metastable vacuum is at the Planck density, the number of trials required is expected to be not too large; but if it is 
at the GUT density which turns out to be many orders of magnitude lower than the Planck density, then the number 
of trials becomes truly formidable (P - lO" 10 ' 8 ) |§. Thus, Farhi, Guth, and Guven H§ guess that it is unlikely 
that the human race will ever succeed in making such a universe in the lab at the GUT scale. Gott applying the 
Copernican principle to estimate our future prospects, would come to similar conclusions. However, if our universe 
is open, it has an infinite number of galaxies, and it would likely have some super-civilizations powerful enough to 
succeed at such a creation event, or at least have so many super-civilizations (an infinite number) that even if they 
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each tried only a few times, then some of them (again an infinite number) would succeed. In fact, if the probability 
for a civilization to form on a habitable planet like the Earth and eventually succeed at creating a universe in the lab 
is some finite number greater than zero (even if it is very low), then our universe (if it is an open bubble universe) 
should spawn an infinite number of such baby universes. 

This notion has caused Harrison |}3| to speculate that our universe was created in this way in the lab by some 
super-civilization in a previous universe. He noted correctly that if super-civilizations in a universe can create many 
baby universes, then baby universes created in this way should greatly outnumber the parent universes, and that 
you (being not special) are simply likely to live in one of the many baby universes, because there are so many more 
of them. Here he is using implicitly the formulation of Gott |6j| that according to the Copernican principle, out of 
all the places for intelligent observers to be, there are, by definition, only a few special places and many non-special 
places, and you are simply more likely to be in one of the many non-special places. Thus, if there are many baby 
universes created by intelligent supercivilizations in an infinite open bubble universe, then you are likely to live in 
a baby universe created in this way. Harrison uses this idea to explain the strong anthropic principle. The strong 
anthropic principle as advanced by Carter p2| says that the laws of physics, in our universe at least, must be such 
as to allow the development of the intelligent life. Why? Because we are here. It is just a self-consistency argument. 
This might lead some to believe, particularly with inflationary cosmologies that are capable of producing an infinite 
number of bubble universes, that these different universes might develop with many different laws of physics, given 
a complicated, many-dimensional inflationary potential with many different minima, and many different low energy 
laws of physics. If some of these did not allow the development of intelligent life and some of these did, well, which 
type of universe would you expect to find yourself in? — one that allowed intelligent observers, of course. (By the 
same argument, you are not surprised to find yourself on a habitable planet — Earth — although such habitable 
planets may well be outnumbered by uninhabitable ones — Mercury, Venus, Pluto, etc.) Thus, there may be many 
more universes that have laws of physics that do not allow intelligent life — you just would not find yourself living 
there. It has been noticed that there are various coincidences in the physical constants — like the numerical value of 
the fine structure constant, or the ratio of the electron to proton mass, or the energy levels in the carbon nucleus — 
which, if they were very different, would make intelligent life either impossible, or much less likely. If we observe such a 
coincidence, according to Carter [p2[ , it simply means that if it were otherwise, we would not be here. Harrison |33| has 
noted that if intelligent civilizations made baby universes they might well, by intelligent choice, make universes that 
purposely had such coincidences in them in order to foster the development of intelligent life in the baby universes 
they created. If that were the case, then the majority of universes would have laws of physics conducive to the 
formation of intelligent life. In this case, the reason that we observe such coincidences is that a previous intelligent 
civilization made them that way. One might even speculate in this scenario that if they were smart enough, they 
could have left us a message of sorts in these dimensionless numbers (a theme that resonates, by the way, with part 
of Carl Sagan's thesis in Contact). However, it is unclear whether any super-civilization would be able to control the 
laws of physics in the universes they created. All, they might reasonably be able to do would be to drive the baby 
universe up into a particular metastable vacuum fl32| . But then, such a metastable vacuum inflates in the knob, and 
an infinite number of bubble universes form later, with perhaps many different laws of physics depending on how they 
tunnel away from the metastable vacuum and which of the many potential minima they roll down into. Controlling 
these phase transitions would seem difficult. Thus, it would seem difficult for the super-civilization that made the 
metastable state that later gave rise to our universe to have been able to manipulate the physical constants in our 
universe. Harrison's model could occur in many generations, making it likely that we were produced as great, great, 
great grandchildren universes from a sequence of intelligent civilizations. Harrison |?3| was able to explain all the 
universes by this mechanism except for the first one! For that, he had to rely on natural mechanisms. This seems to 
be an unfortunate gap. In our scenario, suppose that "first" universe simply turned out to be one of the infinite ones 
formed later by intelligent civilizations. Then the Universe — note capital U — would be multiply connected, and 
would have a region of CTCs; all of the individual universes would owe their birth to some intelligent civilization in 
particular in this picture. 

All this may overestimate the importance of intelligent civilizations. It may be that bubbles of inflating metastable 
vacuum are simply produced at late times in any Big Bang cosmology by natural processes, and that baby universes 
produced by natural processes may vastly outnumber those produced by intelligent civilizations. Such a mechanism 
has been considered by Frolov, Markov, and Mukhanov |5j| . They considered the hypothesis that spacetime curvature 
is limited by quantum mechanics and that as this limit is approached, the curvature approaches that of de Sitter space. 
Then, as any black hole collapses, the curvature increases as the singularity is approached; but before getting there it 
will convert into a collapsing de Sitter solution. This can be done in detail in the following way. Inside the horizon, but 
outside the collapsing star the geometry becomes Schwarzschild which is a radially collapsing but stretching cylinder. 
This can be matched onto a radially collapsing and radially shrinking cylinder in de Sitter space as described by the 
metric in Eq. ( f72[ ) with the time r being negative and the coordinate I being unbounded rather than periodic. Both 
surfaces are cylinders with identical intrinsic curvature, but with different extrinsic curvature. This mismatch is cured 
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by introducing a shell of matter which converts the stretching of the Schwarzschild cylinder to collapsing as well which 
then matches onto the collapsing de Sitter solution. This phase transition may occur in segments which then merge 
as noted by Barabes and Frolov [^6|j5^]. The de Sitter solution then bounces and becomes an expanding de Sitter 
solution which can in turn spawn an infinite number of open bubble universes. This all happens behind the event 
horizon of the black hole. Within the de Sitter phase, one finds a Cauchy horizon like the interior Cauchy horizon 
of the Reisner-Nordstrom solution, but this inner Cauchy horizon is not unstable because the curvature is bounded 
by the de Sitter value so the curvature is not allowed to blow up on the inner horizon. (This is an argument that 
one could also rely on to produce self-consistent multiply connected de Sitter phases with CTCs — if needed.) This 
model thus produces, inside the black hole, to the future, and behind the event horizon, an expanding de Sitter phase 
that has a beginning, just like Vilenkin's tunneling universe. If one of those bubble universes simply turns out to be 
the original one in which the black hole formed, then the solution is multiply connected with a region of CTCs. This 
would make every black hole produce an infinite number of universes. This would be the dominant mechanism for 
making new bubble universes, since the number of black holes in our universe would appear to greatly outnumber the 
number of baby universes ever produced by intelligent civilizations, since the tunneling probability for that process 
to succeed is exceedingly small. 

Smolin 1 34 3f| has proposed that this type of mechanism works and furthermore that the laws of physics (in the 
bubble universes) are like those in our own but with small variations. Then, there would be a Darwinian evolution of 
universes. Universes that produced many black holes would have more children that would inherit their characteristics 
— with some small variations. Soon, most universes would have laws of physics that were fine-tuned to produce the 
maximum number of black holes. Smolin |34||35| points out that this theory is testable, since we can calculate whether 
small changes in the physical constants would decrease the number of black holes formed. In this picture we should 
be near a global maximum in the black hole production rate. One problem is that the laws of physics that maximize 
the number of black holes and those that simply maximize the number of main sequence stars may be rather similar, 
and the laws that maximize the number of main sequence stars might well simply maximize the number of intelligent 
observers, and the anthropic principle alone would suggest a preference for us observing such laws, even if no baby 
univ erses were created in black holes. Another possible problem with this model, pointed out by Rothman and Ellis 
[186|, is that if the density fluctuations in the early universe had been higher in amplitude, this would form many tiny 
primordial black holes (presumably more black holes per comoving volume than in our universe), so, we well might 
wonder why the density fluctuations in our universe were so small. One way out might be that tiny black holes do 
not form any baby universes, but this seems a bit forced since the de Sitter neck formed can be as small as the Planck 
scale or GUT scale and it would seem that even primordial black holes could be large enough to produce an infinite 
number of open bubble universes. 

Another possibility is the recycling universe of Garriga and Vilenkin [ p6[ . In this model there is a metastable 
vacuum with cosmological constant Aj., and a true lowest energy vacuum with a cosmological constant A2. Ai is at 
the GUT or Planck energy scale, while A2 is taken to be the present value of A (as might be the case in a flat-A 
model). As long as A2 > 0, then Garriga and Vilenkin assert that there is a finite (but small) probability per unit 
four volume that the A2 state could tunnel to form a bubble of Ai state, which could therefore inflate, decaying into 
bubbles of A2 vacuum, which could recycle forming Ai bubbles, and so forth. They point out that depending on the 
coordinate system, a bubble of A2 forming inside a Ai universe could also be seen as a Ai bubble forming inside of a 
A2 universe. Take two de Sitter spaces, one with Ai and one with A2, and cut each along a vertical slice (W = Wq) 
in the embedding space. They can then be joined along an appropriate hyperbola of one sheet representing a bubble 
wall, with the A 2 universe lying to the W < Wq side and the Ai universe lying to the X > Wo side. Slicing along 
hyperplanes with V + W = constant gives a steady-state coordinate system for a Ai universe in which a bubble of A2 
vacuum appears. Slicing along hyperplanes with V — W = constant, however, gives a steady-state coordinate system 
for a A2 universe in which a bubble of Ai appears. So, one can find a steady-state coordinate system in which there 
is a Ai universe, with bubbles of A2 inside it, and bubbles of Ai inside these A2 bubbles, and so forth. If the roll 
down is slow, within the A2 bubble as it forms, as in Gott's open bubble universe |^{|, then it will have at least 67 
e- folds of inflation with A ~ Ai before it falls off the plateau into the absolute minimum at A2, and this will be an 
acceptable Big Bang model which will have the usual Big Bang properties except that it will eventually be dominated 
by a lambda term A2. Being bubble universes, they will all be open with negative curvature as in Gott's model |8^] 
but they will be asymptotically open de Sitter models at late times with a(t) = ro sinh(t/ro) and A = A2. Garriga 
and Vilenkin |36f| wondered whether such a recycling model could be geodesically complete toward the past. Such a 
outcome, they pointed out, would violate no known theorems and should be investigated. They hoped to find such a 
geodesically-complete-to-the-past model so that it could be eternal without a need for a beginning. However, in the 
special case, where Ai = A2, one can show that the recycling steady state solution becomes a simple single de Sitter 
space geometry with Ai and the usual steady-state coordinate system in a single de Sitter space is not geodesically 
complete to the past. 

Now take this recycling model where it turns out that one of the Ai bubbles formed inside an A2 bubble inside 
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a Ai region is, in fact, the Ai region that one started out with. In this case, we would have a multiply connected 
model such as we are proposing which would include a region of CTCs (Fig. g). (If Ai = A2, this model is just the 
multiply connected de Sitter space we have considered.) If our multiply connected model was geodesically complete 
to the past, so would the covering space (a simply connected Garriga-Vilenkin model) be. If our multiply connected 
model was geodesically incomplete to the past, so would the covering space (a simply connected Garriga-Vilenkin 
model) be also. In our model, there would be a strong self-consistency reason for pure retarded potential, whereas 
in the Garriga-Vilenkin recycling model, there would be no such strong reason for it. With pure retarded potentials 
throughout, the issue of whether the spacetime was geodesically complete to the past is less compelling, as we have 
argued above, and our model, having a periodic boundary condition, would not need further boundary conditions, 
unlike a simply connected recycling model that was geodesically incomplete to the past. 

Thus, there are a number of models in which baby universes are created which can be converted into models in 
which the Universe creates itself, if one of those created baby universes turns out to be the original universe that 
one started with. Since these models are all ones in which there are an infinite number of baby universes created, 
this multiply connected outcome must occur unless the probability for a particular multiple connectivity to exist is 
exactly zero. In other words, it should occur, unless it is forbidden by the laws of physics. Given quantum mechanics, 
it would seem that such multiple connectivities would not be absolutely forbidden, particularly in the Planck foam 
era. 

We should note here that, in principle, there might even be solutions that are simply connected in which there was 
an early region of CTCs bounded to the future by a Cauchy horizon followed by an inflationary region giving rise to 
an infinite number of bubble universes. The models considered so far have all obeyed the weak energy condition, and 
these models have all been multiply connected; in other words, they have a genus of 1, like a donut, since one of the 
later baby universes is connected with the original one. Consider an asymptotically flat spacetime with two connected 
wormhole mouths that are widely separated. The existence of the wormhole connection increases the genus by one. 
Instead of a flat plane, it becomes a flat plane with a handle. To do this, the wormhole solution must violate the weak 
energy condition |Q. It must have some negative energy density material, for it is a diverging lens (converging light 
rays entering one wormhole mouth, diverge upon exiting the other mouth). For a compact two dimensional surface, 
the integrated Gaussian curvature over the surface divided by An is equal to 1 minus the genus. Thus, the integrated 
Gaussian curvature over a sphere (genus=0) is An, while the integrated Gaussian curvature over a donut (genus=l) 
is zero, and the integrated Gaussian curvature over a figure 8 pretzel (genus=2) is —An. Negative curvature is added 
each time the genus is increased. Conversely, positive curvature can be added to reduce the genus by 1. When a donut 
is cut, so that it resembles a letter "C" , the ends of the letter "C" are sealed with positive curvature (two spherical 
hemispherical caps would do the job, for example). Our solutions are already multiply connected, so they might in 
principle be made simply connected by the addition of some extra positive mass density, without violating the weak 
energy condition. An example of this is seen by comparing Grant space |L8| with Gott's two-string spacetime [||. 
Grant space is multiply connected, has T a b = everywhere, and includes CTCs. It can be pictured as a cylinder. 
Gott's two-string spacetime is simply connected, but is identical to Grant space at large distances from the strings. 
It also contains CTCs. It can be pictured as a cardboard cylinder that has been stepped on and then stapled shut 
at one end, like an envelope. There are two corners at the closed end, representing the two strings, but the cylinder 
continues outward forever toward its open end (so it is like a test tube, a cylinder closed on one end). The two strings 
provide positive energy density (i.e. they do not violate the weak energy condition). CTCs that wrap around the two 
strings far out in the cylinder (which is identical to a part of Grant space; see Laurence [187Q can be shrunken to 
points by slipping them through the strings — but they become spacelike curves during this process. Thus, Gott space 
represents how a multiply connected spacetime with CTCs (Grant space) can be converted into a simply connected 
spacetime with CTCs by adding to the solution material that obeys the weak energy condition. A similar thing might 
in principle be possible with these cosmological models. Since our multiply connected versions already obey the weak 
energy condition, so would the associated simply connected versions. 



XI. CONCLUSIONS 



The question of first-cause has been a troubling one for cosmology. Often, this has been solved by postulating a 
universe that has existed forever in the past. Big Bang models supposed that the first-cause was a singularity, but 
questions about its almost, but not quite, uniformity remained. Besides, the Big Bang singularity just indicated a 
breakdown of classical general relativity, and with a proper theory-of-everything, one could perhaps push through 
to earlier times. Inflation has solved some of these problems, but Borde and Vilenkin have shown that if the initial 
inflationary state is metastable, then it must have had a finite beginning also. Ultimately, the problem seems to be 
how to create something out of nothing. 
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So far, the best attempt at this has been Vilenkin's tunneling from nothing model and the similar Hartle-Hawking 
no-boundary proposal. Unfortunately, tunneling is, as the name suggests, usually a process that involves tunneling 
from one classical state to another, thus, with the Wheeler-DeWitt potential and "energy" E = that Hartle and 
Hawking adopted, the Universe, we argue, should really start not as nothing but as an S 3 universe of radius zero — 
a point. A point is as close to nothing as one can get, but it is not nothing. Also, how could a point include the laws 
of physics? In quantum cosmology, the wave function of the Universe is treated as the solution of a Schrodinger-likc 
equation (the Wheeler-DeWitt equation), where the three-sphere S 3 radius a is the abscissa and there is a potential 
U(a) with a metastable minimum at U(a = 0) = 0, and a barrier with U(a) > for < a < do, and U(a) < for 
a > oq. Thus, the evolution can be seen as a particle, representing the universe, starting as a point, a = 0, at the 
bottom of the metastable potential well, with E = 0. Then it tunnels through the barrier and emerges at a = ao 
with E — 0, whereupon it becomes a classically inflating de Sitter solution. It can then decay via the formation 
of open single bubble universes |89||)8]]. The problem with this model is that it ignores the "zero-point-energy". If 
there is a conformal scalar field </>, then the "energy" levels should be E n = n + h. Even for n = there is a "zero- 
point-energy". The potential makes the system behave like a harmonic oscillator in the potential well near a = 0. 
A harmonic oscillator cannot sit at the bottom of the potential well — the uncertainty principle would not allow 
it. There must be some zero-point-energy and the particle must have some momentum, as it oscillates within the 
potential well when the held <f> is included. Thus, when the "zero-point-energy" is considered, we see that the initial 
state is not a point but a tiny oscillating (0 < a < a±) Big Bang universe, that oscillates between Big Bangs and Big 
Crunches (though the singularities at the Big Bangs and Big Crunches might be smeared by quantum effects). This 
is the initial classical state from which the tunneling occurs. It is metastable, so this oscillating universe could not 
have existed forever: after a finite half-life, it is likely to decay. It reaches maximum radius a\, and then tunnels to 
a classical de Sitter state at minimum radius 02 where a-i < a . The original oscillating universe could have formed 
by a similar tunneling process from a contracting de Sitter phase, but such a phase would have been much more 
likely to have simply classically bounced to an expanding de Sitter phase instead of tunneling into the oscillating 
metastable state at the origin. In this case, if one found oneself in an expanding de Sitter phase, it would be much 
more likely that it was the result of classical bounce from a contracting de Sitter phase, rather than the result of a 
contracting de Sitter phase that had tunneled to an oscillating phase and then back out to an expanding de Sitter 
phase. Besides, a contracting de Sitter phase would be destroyed by the formation of bubbles which would percolate 
before the minimum radius was ever reached. 

In this paper, we consider instead the notion that the Universe did not arise out of nothing, but rather created 
itself. One of the remarkable properties of the theory of general relativity is that in principle it allows solutions with 
CTCs. Why not apply this to the problem of the first-cause? Usually the beginning of the Universe is viewed like the 
south pole. Asking what is before that is like asking what is south of the south pole, it is said. But as we have seen, 
there remain unresolved problems with this model. If instead there were a region of CTCs in the early universe, then 
asking what was the earliest point in the Universe would be like asking what is the easternmost point on the Earth. 
There is no easternmost point — you can continue going east around and around the Earth. Every point has points 
that are to the east of it. If the Universe contained an early region of CTCs, there would be no first-cause. Every 
event would have events to its past. And yet the Universe would not have existed eternally in the past (see Fig. ^). 
Thus, one of the most remarkable properties of general relativity — the ability in principle to allow CTCs — would 
be called upon to solve one of the most perplexing problems in cosmology. Such an early region of CTCs could well 
be over by now, being bounded to the future by a Cauchy horizon. We construct some examples to show that vacuum 
states can be found such that the renormalized energy-momentum tensor does not blow up as one approaches the 
Cauchy horizon. For such a model to work the Universe has to reproduce at some later time the same conditions that 
obtained at an earlier time. Inflation is particularly useful in this regard, for starting with a tiny piece of inflating 
state, at later times a huge volume of inflating state is produced, little pieces of which look just like the one we started 
with. Many inflationary models allow creation of baby inflationary universes inside black holes, either by tunneling 
across the Einstein-Rosen bridge, or by formation as one approaches the singularity. If one of these baby universes 
simply turns out to be the universe we started with, then a multiply connected model with early CTCs bounded by 
a Cauchy horizon is produced. Since any closed null geodesies generating the Cauchy horizon must circulate through 
the optically thick region of the hot Big Bang phase of the universe after the inflation has stopped, the renormalized 
energy-momentum tensor should not blow up as the Cauchy horizon is approached. 

As a particularly simple example we consider a multiply connected de Sitter solution where events Ej are topologi- 
cally identified with events E^ that lie inside these future light cones via a boost transformation. If the boost b = 2tt, 
we show that we can find a Rindler-type vacuum where the renormalized energy-momentum tensor does not blow up 
as the Cauchy horizon is approached but rather produces a cosmological constant throughout the spacetime which 
self-consistently solves Einstein's equations for this geometry. Thus, it is possible to find self-consistent solutions. 
When analyzing classical fields in this model, the only self-consistent solution without a blow up as the Cauchy hori- 
zon is approached occurs when there is a pure retarded potential in the causally connected region of the model. Thus, 
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the multiply connected nature of this model and the possibility of waves running into themselves, ensure the creation 
of an arrow of time in this model. This is a remarkable property of this model. Interestingly, this model, although 
having no earliest event and having some timelike geodesies that are infinitely extendible to the past, is nevertheless 
geodesically incomplete to the past. This is not a property we should have thought desirable, but since pure retarded 
potentials are established automatically in this model, there are no waves propagating to the past and so there may 
be no problem in physics with this, since there are never any waves that run off the edge of the spacetime. The region 
of CTCs has a finite four- volume equal to A-wbr^ /3 and should be in a pure vacuum state containing no real particles 
or Hawking radiation and no bubbles. After the Cauchy horizon for a certain amount of proper time (depending on 
the bubble formation probability per four volume fg) no bubbles (or real particles) form, but eventually this model 
expands to infinite volume, creating an infinite number of open bubble universes, which do not percolate. At late 
times in the de Sitter phase a particle detector would find the usual Hawking radiation just as in the usual vacuum 
for de Sitter space. 

There are a number of problems to be solved in this model. The chronology projection conjecture proposes that 
the laws of physics conspire so as to prevent the formation of CTCs. This conjecture was motivated by Hiscock and 
Konkowski's result that the energy-momentum tensor of the adapted Minkowski vacuum in Misner space diverges as 
the Cauchy horizon is approached. But as we have shown ]30[] , the adapted Rindler vacuum for Misner space has 
(Tat) rcn = throughout the space if b — 2n; thus, this is a self-consistent vacuum for this spacetime since it solves 
Einstein's equations for this geometry. It's true that (T a f,) ren remains formally ill-defined on the Cauchy horizon itself 
[£ = in Eq. ( |32] ) with b = 2tt], a set of measure zero. But it is not clear that this creates a problem for physics, 
since continuity might require that this formally ill-defined quantity be defined to be zero on this set of measure zero 
as well, since it is zero everywhere else. In fact, a treatment in the Euclidean section shows this is the case, for in the 
Euclidean section, if b = 2tt, (T a b) len — everywhere, including at £ = 0. Other counter-examples to the chronology 
protection conjecture have also been found, as discussed in section V. Hawking himself has also admitted that the 
back-reaction of vacuum polarization does not enforce the chronology protection conjecture. 

One of the remarkable properties of general relativity is that it allows, in principle, the formation of event horizons. 
This appears to be realized in the case of black holes. Just as black hole theory introduced singularities at the end, 
standard Big Bang cosmology introduced singularities at the beginning of the universe. Now, with inflation, we see 
that event horizons should exist in the early universe as well [|8{| . Inflationary ideas prompt the suggestion that baby 
universes may be born. If one of the baby universes simply turns out to be the one we started with, then we get a 
model with an epoch of CTCs that is over by now, bounded toward the future by a Cauchy horizon. We have argued 
that the divergence of the energy-momentum tensor as one approaches the Cauchy horizon does not necessarily occur, 
particularly when the Cauchy horizon crosses through a hot Big Bang phase where absorption occurs. 

If the energy-momentum tensor does not diverge as the Cauchy horizon is approached, other problems must still 
be tackled. The classical instability of a Cauchy horizon to the future (a future chronology horizon) in a spacetime 
with CTCs is one. But this problem is solved in a world with retarded potentials for a Cauchy horizon that occurs 
to our past (a past chronology horizon) and which ends an epoch of CTCs. It thus seems easier to have a Cauchy 
horizon in the early universe. At the microscopic level, quantum mechanics appears to allow acausal behavior. Indeed 
the creation and annihilation of a virtual positron-electron pair can be viewed as creation of a small closed loop, 
where the electron traveling backward in time to complete the loop appears as a positron. So, why should the laws 
of physics forbid time travel globally? Indeed one of the most remarkable properties of the laws of physics is that 
although they are time (CPT) symmetric, the solutions we observe have an arrow of time and retarded potentials. 
Without this feature of the solutions, acausal behavior would be seen all the time. Interestingly, in our model, the 
multiply connected nature of the spacetime geometry forces an arrow of time and retarded potentials. Thus, it is the 
very presence of the initial region of CTCs that produces the strong causality that we observe later on. This is a very 
interesting and unexpected property. An entropy arrow of time is automatically produced as well, with the region of 
CTCs in the simplest models sitting automatically in a co ld vacuum state, with the universe becoming heated after 
the Cauchy horizon. Recently, Cassidy and Hawking [142] have proposed yet another supposed difficulty for CTCs, 
in that the formally defined entropy appears to diverge to negative infinity as the Cauchy horizon is approached. Yet, 
in the early universe this may turn out to be an advantage, since to produce the ordinary entropy arrow of time we 
observe in the universe today, we must necessarily have some kind of natural low-entropy boundary condition in the 
early universe |5^,|5^]. This could occur on the Cauchy horizon that ends the period of CTCs. 

New objections to spacetimes with CTCs can continue to surface, as old problems are put to rest, so it might seem 
that disproving the chronology protection conjecture would be a tall order. But, proving that there are no exceptions 
to the chronology protection conjecture, ever, would seem an equally daunting task. This is particularly true since 
we currently do not have either a theory of quantum gravity or a theory-of-everything. 

Perhaps the most obvious problem with the model we have proposed is that the simplest solutions we have obtained 
so far are not geodesically complete to the past. But we may need no boundary condition since we have a periodic 
boundary condition instead. This thus may not be a problem in physics if retarded potentials are the only ones 
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allowed. Alternatively, as Garriga and Vilenkin have indicated, it would violate no known theorems for some type 
of recycling universe (making bubble universes within bubble universes ad infinitum) to exist that was geodesically 
complete to the past. If such solutions exist, it might be possible to find a solution in which there was an early epoch 
of CTCs that would be geodesically complete to the past as well by simply identifying an earlier bubble with a later 
one. 

Thus, a number of important questions remain, and we would not minimize them. The models presented here, 
however, do have some interesting and attractive properties, suggesting that this type of model should be investigated 
further, and that we should ask the question: 

Do the laws of physics prevent the Universe from being its own mother? 
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FIG. 1. A self-creating Universe scenario. Four inflating baby universes are shown — A, B, C, and D — from left to right. 
Universe A and D have not created any baby universes so far. Universe C has created universe D. Universe B has created three 
universes: universe A, universe C and itself. The toroidal — shaped region at the bottom is a region of CTCs (closed timelike 
curves). The region is bounded to the future by a Cauchy horizon, after which, there are no CTCs. Universes A, C, and D, for 
example, are formed after the Cauchy horizon when the epoch of CTCs is already over. 
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FIG. 2. The potential function in the Wheeler-DeWitt equation in the minisuperspace model. The horizontal axis is the 
scale factor of the universe. If the conformally coupled scalar field is in the ground state, it has a "zero-point-energy" 1/2. If 
this "zero-point-energy" is considered, the quantum behavior of the universe is like a particle of unit mass with total energy 
1/2 moving in the potential {7(a). Regions < a < ai and a > a?, are classically allowed, region a 1 < a < ai is classically 
forbidden. The left dark disk is a tiny radiation-dominated closed oscillating universe, which oscillates between Big Bangs 
and Big Crunches. The smoothness of the potential at a = may indicate that any Big Bang and Big Crunch singularities 
are removed by quantum theory. This tiny oscillating universe has a small but non-zero probability to tunnel through the 
barrier out to become a de Sitter-type inflating universe, which is represented by the dark disk on the right. The circle inside 
the barrier is a Euclidean bouncing space. If the "zero-point-energy" 1/2 were neglected (as Hartle and Hawking did), the 
left classically allowed region would shrink to a point. The grey disk represents a contracting and re-expanding de Sitter 
universe. If the "zero-point-energy" is neglected, the Universe could start out at the metastable minimum as a point with 
a — 0, tunneling through the barrier out to become a de Sitter universe. In this paper we argue that we have no reason to 
neglect the "zero-point-energy" so that it is the tiny oscillating universe initial state that applies. 

FIG. 3. (a) The solution of the Euclidean Einstein's equations representing the tunneling regime (open circle) in Fig. ^| This 
is a solution to the Euclidean Einstein's equations with a positive cosmological constant and a conformally coupled scalar field 
in its ground state. This is a Euclidean space bouncing between the state with maximum radius a% and the state with (non-zero) 
minimum radius 01. One "copy" of this Euclidean bouncing solution is shown in this diagram, which has two boundaries with 
minimum radius 01. (b) This is the case when the "zero-point-energy" of the conformally coupled scalar field is neglected, as 
Hartle and Hawking did. In this case the minimum radius is zero, and thus one copy of the bouncing Euclidean solution is a 
four-sphere. This four-sphere has no-boundary, which is the basis of Hawking's quantum cosmology. But we argue that since 
the "zero-point-energy" of the conformally coupled scalar field cannot be neglected, the true solution should be that given by 
diagram 3a, which does not enforce Hartle and Hawking's no-boundary proposal. 

FIG. 4. Penrose diagrams of our multiply connected de Sitter space mapped onto its universal covering space (de Sitter space). 
Under a boost transformation, points with the same symbols (squares, disks, triangles, or double-triangles) are identified. Our 
multiply connected de Sitter space is divided into four regions TZ, C, and V, which are separated by Cauchy horizons CTL. 
The shaded regions represent fundamental cells of the multiply connected de Sitter space. (Fig. ^a and Fig. ^b represent two 
different choices of the fundamental cells, but they are equivalent.) The fundamental cells TZ and C have a finite four-volume, 
whereas the fundamental cells T and V (which extend infinitely to the future and the past, respectively) have an infinite 
four-volume. In Fig. ^a the left and right boundaries of T are identified, likewise for V; the upper and lower boundaries of TZ 
are identified, likewise for £. In Fig. ^b region T + TZ is partially bounded by two null surfaces, the lower one is the future light 
cone of an event E, and the upper one is the future future light cone of an event E' which is identified with E under the action 
of a boost. These two future light cones are identified creating a periodic boundary condition for the causally connected region 
J- + TZ. 7Z and T are separated by a Cauchy horizon CTL. Self-consistency (non-divergence of T ab T ab as CTL is approached) 
requires retarded potentials in 1Z and T . Region V + C is partially bounded by the past light cone of an event F and the past 
light cone of of an event F' which is identified with F under the action of a boost. These two surfaces are identified creating 
a periodic boundary condition for V + L, where self-consistency as CTL separating V from C is approached requires advanced 
potentials. 

FIG. 5. With our adapted conformal Rindler vacuum, our multiply connected de Sitter space is cold (with zero temperature) 
in 1Z and C, but hot (with the Gibbons-Hawking temperature) in T and V ■ The arrows indicate the direction of increasing 
entropy. 

FIG. 6. Self-consistency near the Cauchy horizons in a spacetime with CTCs naturally gives rise to an arrow of time. Grey 
thick lines represent light cones of electromagnetic waves or photons emitted from event E. (a) This diagram shows that in T 
the retarded potential is self-consistent. The "collision" of an electromagnetic wave with its images cannot destroy the Cauchy 
horizon, since the proper time from the "collision" (event p) to the origin is always bigger than the proper time from E to the 
origin. Likewise the advanced potential in region V would not destroy the Cauchy horizon, (b) This diagram shows that a 
retarded potential in TZ and an advanced potential in C (or vice versa) are self-consistent. But the potentials in 1Z and C cannot 
be both retarded or both advanced, otherwise the "collision" of two waves from TZ and £ respectively will destroy the Cauchy 
horizon, (c) This diagram shows that the advanced potential in T is not self-consistent, since the collision of an electromagnetic 
wave with its images will destroy the Cauchy horizon, (as n — » ±00, the collision event p approaches the Cauchy horizon.) (d) 
This diagram shows that a part-advanced-and-part-retarded potential in TZ (or £) is also not self-consistent, as T ab T a t would 
also diverge as the Cauchy horizon is approached. 
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FIG. 7. A schematic Penrose diagram of a self-creating Universe based on the baby universe model of Farhi, Guth, and 
Guven. We identify M\N\ with M2-/V2, to obtain a model of the Universe creating itself. (M1N1 is the future light cone of 
event Ni.) In this model the closed null curves generating the Cauchy horizon (CTi.) pass through a hot Big Bang region, where 
the dense absorber can make the Cauchy horizon stable against vacuum polarization effects. The metastable de Sitter phase 
is shown in grey. It decays along a hyperboloid H 3 near the bottom of the figure to form a single open bubble universe with 
a hot Big Bang phase and an epoch of recombination which is also shown. After recombination a super-civilization creates, at 
the right, an expanding bubble of de Sitter metastable vacuum. This reaches a point of maximum expansion at which point it 
tunnels to a doorknob-shaped configuration. The tunneling epoch is shown by the dashed line: just below the dashed line is 
how the spacetime appears just before tunneling, and just above the dashed line is how the spacetime appears just after the 
tunneling. Just after the tunneling, the geometry (just above the dashed line) from left to right goes from infinite radius (where 
future null infinity I + meets the dashed line) to a minimum radius r = 2M at the neck in the Einstein-Rosen bridge (where the 
inside and outside black hole event horizons meet just below the word "Black Hole") then to a radius r > 2AI at the surface 
of the de Sitter bubble, reaching a maximum radius at the equator of the de Sitter bubble "knob" and finally decreasing to 
r = at the center of the bubble at the extreme right. The de Sitter bubble expands forever. M2 is at t — 00. To the left of 
M2 is another open bubble universe forming out of the metastable de Sitter vacuum. It is diamond-shaped — the bottom two 
lines representing the expanding bubble wall and the top two lines representing future null infinity for that bubble. Within 
this bubble the de Sitter vacuum decays to a hot Big Bang phase along a hyperboloid H 3 shown as a curved line crossing the 
diamond. Another open bubble universe forms to the right of Mi. Recall, Mi = M2. These two bubble universes both form 
after the Cauchy horizon CTi as do an infinite number of others. The black hole singularity is shown, as well as the fact that 
the black hole evaporates. N1N2 is a CTC; CTCs occur on the N1N2 side of the Cauchy horizon CTi. After CTi, there are no 
CTCs. 



FIG. 8. A self-creating Universe model based on Garriga and Vilenkin's recycling Universe. In a region of cosmological 
constant Ai , a bubble B of cosmological constant A2 is formed by tunneling at the epoch BBi . The expanding bubble wall is 
represented by BB2- At a later time, within bubble B a bubble A forms at epoch AA\ by tunneling. The expanding bubble 
wall is shown by AA2. Inside bubble A the cosmological constant is Ai. In the limit where Ai = A2 we can plot this in a single 
de Sitter space. Now we identify the two hypersurfaces denoted by A1AA2 to obtain a model of the Universe creating itself. 
The Cauchy horizon bounding the region of CTCs is indicated by CTi. After CTi there are no CTCs. If Ai = A2, this reduces 
to our multiply connected de Sitter model J- +1Z shown in Fig. kp. 
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